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Abstract 

The main goal of this work consists in showing that the analytic solu- 
tions for a class of characteristic problems for the Einstein vacuum equa- 
tions have an existence region much larger than the one provided by the 
Cauchy-Kowalevski theorem due to the intrinsic hyperbolicity of the Ein- 
stein equations. To prove this result we first describe a geometric way of 
writing the vacuum Einstein equations for the characteristic problems we 
are considering, in a gauge characterized by the introduction of a double 
null cone foliation of the spacetime. Then we prove that the existence 
region for the analytic solutions can be extended to a larger region which 
depends only on the validity of the apriori estimates for the Weyl equa- 
tions, associated to the "Bel- Robinson norms" . In particular if the initial 
data are sufficiently small we show that the analytic solution is global. 
Before showing how to extend the existence region we describe the same 
result in the case of the Burger equation, which, even if much simpler, nev- 
ertheless requires analogous logical steps required for the general proof. 
Due to length of this work, in this paper we mainly concentrate on the 
definition of the gauge we use and on writing in a "geometric" way the 
Einstein equations, then we show how the Cauchy-Kowalevski theorem 
is adapted to the characteristic problem for the Einstein equations and 
we describe how the existence region can be extended in the case of the 
Burger equation. Finally we describe the structure of the extension proof 
in the case of the Einstein equations. The technical parts of this last result 
is the content of a second paper. 
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1 Introduction 



In this paper we prove a result about the existence region of the analytic so- 
lution of a class of characteristic problems, namely those whose "initial data" 
are given on a null hypersurface consisting of the union of a truncated outgoing 
null cone and of a truncated incoming cone intersecting the previous one along 
a surface diffcomorphic to S 2 . This class of characteristic problems has been 
studied by different authors, for instance H.Mullcr Zum Hagen, |Muj . H.Mullcr 
Zum Hagen and H.J.Seifert, |Mu-Sej . Christudoulou and Muller Zum Hagen, 
|Ch-MZj . Dossa in a series of papers, see |Do] and references therein, but, in 
particular, we recall the anticipating work by A.Rendall, [Ren], where a thor- 
ough examination has been done to show how to obtain initial data satisfying 
the costraint equations and the harmonic gauge conditions and, subsequently, 
a way of obtaining a local existence result is presented. Recently following, but 
largely improving the A.Rendall result, we suggest to the reader attention the 
paper by Y. Choquet-Bruhat, P.Y.Chrusciel, J.M. Martin-Garcia, "The Cauchy 
problem on a characteristic cone for the Einstein equations in arbitrary dimen- 
sions" , to appear on Annales Henry Poincare . In that paper the authors prove 
a local existence result for the characteristic problem with initial data on a null 
cone, using again the harmonic gauge and proving in a very detailed way how 
the initial data constraints have to be satisfied and how, relying on the Dossa 
results, the local existence result can be proved. Moreover the nature of the 
characteristic problem, initial data on the null cone, adds the extra problems of 
the "tip of the cone" they solve completely. 

In the present paper our goal is to show that the real analytic solutions of the 
class of characteristic problems we are considering have, due to the hyperbolicity 
of the Einstein equations, a larger existence region than the one we can obtain 
by the application of the Cauchy-Kowalevski theorem. More precisely we prove 
that the extension of the analyticity existence region depends only on a finite 
number of derivatives of the initial data, namely on some appropriate Sobolcv 
norms; moreover if we assume the initial data small in these norms we can 
prove the global existence of the analytic solutions. Some analogous results 
have been obtained in the past by S.Alinhac, G.Metivier, where they proved 
the propagation of the analiticity for hyperbolic systems of p.d.e., see [A-M] and 
references therein; more recent results can be found in S.Spagnolo, "Propagation 
of analyticity for a class of nonlinear hyperbolic equations" , |Spg , and references 
therein. To prove this result we present a strategy analogous to that used 
by S.Klainerman and one of the authors (F.N.) in |K1-Nil] . for the Burger 
equation in a non characteristic problem. Clearly, as the Einstein equations are 
much more complicated, the extension to the present case is significantly more 
difficult. Nevertheless some general aspects are borrowed from that toy model 
and suggests the nature of the different technical problems we have to deal. 
As it will appear clearly in the rest of this introduction, in the bulk of this 
paper and of the following one, preliminary to the existence proof of the analytic 
solutions a detailed examination of the gauge used and how the contraints are 
satisfied is needed. This part, crucial to the development of the existence proof, 
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is, in our opinion, of interest in itself and new in many aspects. To it a greater 
part of the present paper is devoted. Let us summarize the various steps of our 
approach. 

i) As discussed in detail in Subsection 14.11 and stated before the crucial ingre- 
dient to extend the analiticity region is the hyperbolicity of the p.d.e. equations 
we are considering. The region where the analitical solution exists is the region 
where the apriori energy estimates can be proved. Therefore the control of en- 
ergy norms is the key step. To achieve these norms and their apriori estimates 
is trivial in the Burger case, but much more delicate in the Einstein equations, 
even more if we look for energy type estimates valid (for small initial data) ev- 
erywhere. It is, in fact, well known that the standard energy norms associated 
to the vacuum Einstein equation in the harmonic gauge are very difficult to use 
to obtain a global solution even in the noncharacteristic problem, see [Lin-Rod] , 
and no results are at our disposal for the characteristic case. It turns out that an 
efficient strategy to achieve a global existence result is the one introduced first 
by D. Christoudoulou and S. Klainerman in [Ch-Kl] and subsequently modified 
by S. Klainerman and F. Nicolo, see [K1-M2]. In this approach the "energy type" 
norms to bound are those associated to the Bel-Robinson tensor, quadratic in 
the Riemann tensor. The control of these norms is strictly tied to the control 
of the connection coefficient;!]] of the spacetime and the equations which control 
these last quantities are the so called "structure equations", see |Sp| , Vol 2, 
which have the form of elliptic Hodge systems or of transport equations along 
the null directions. 

ii) This strongly suggests the use of a foliation already introduced in |K1-Ni2j , 
Chapter 3, the "double null cone foliation". The main differences are that first, 
the structure equations were used there to obtain good estimates for the var- 
ious L p integral norms while here to show how the Einstein equations can be 
expressed as a subset of these structure equations. More precisely, as in the non 
characteristic case one can foliate the spacetime with a family of spacelike hy- 
persurfaces and write the Einstein equations as a system of first order equations 
for the (riemannian) metric and the second fundamental form adapted to this 
foliation, here we assume the spacetime foliated by a family of outgoing cones 
and incoming (truncated) cones and write the Einstein equations as a set of first 
order equations involving the metric adapted to these cones and the connection 
coefficients, basically, the first derivatives of this metric. The complete detailed 
description of the procedure to write the Einstein equations in the way we are 
sketching here is given in subsections 12.31 and 12.41 

The second difference is that in [K1- M2] a local solution was already assumed 
to exist and proved in the more standard harmonic gauge, here even the local 
existence is proved in the "double null cone foliation gauge" . Therefore in this 
approach we will never use any foliation made by spacelike hypersurfaces, the 
derivatives of the various unknown functions of our equations are always done 
with respect to the angular variables and to the u and u variables, the affine 
parameters of the null geodesies generating the outgoing or incoming cones. 

1 Sometimes called Ricci coefficients. 
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iii) Exactly as in the non characteristic case the choice of the spacelike 
hypersurface foliation is basically equivalent to the choice of a gauge, in the 
present case specifying the "double null cone foliation" is just the choice of the 
gauge and the quantities fl and X which appear in the expression of our metric 
in the adapted coordinates, 

g = -2n 2 (dudu+dudu)+jab(X a du + duj a )(X b du + duj b ) , (1.1) 

play the role of the lapse function and the shift vector. Analogously to them 
they will have to satisfy some differential equations. In subsection 12.31 a very 
extended discussion about this gauge is given. 

iv) As we are dealing with a characteristic problem it is expected that the 
initial data cannot be given in a complete free way, but they have to satisfy 
some constraints. In the case of the Einstein equations the situation is more 
complicated as even in the non characteristic case the initial data cannot be 
given in a free wayH Therefore our initial data have to satisfy two different kind 
of constraints, those due to the nature of the Einstein equations, the analogous 
of the constraints equations for the second fundamental form fey and those 
associated to the gauge choice namely, in the present case, the equations for f2 
and X. In our presentation, more geometric than the one using, for instance, 
the harmonic gauge, the distinction between these two kinds of constraints is 
completely clear and it is natural that in these equations no transverse (to the 
cones) derivatives appear. 

Finally to show that the solutions of our equations are in fact solutions of the 
Einstein equations we have to prove, exactly as in the non characteristic case, 
that the (Einstein) constraint equations once satisfied by the initial data are 
satisfied everywhere. This is proved in subsection 12.41 

v) As we want to prove that our analytic solutions can be extended to the 
whole spacetim^f| we have first to provide a local analitic solution of the char- 
acteristic problem. This is discussed in Section [3] where we adapt the Cauchy- 
Kowalevski theorem to the characteristic problem following G.F.D.Duff, |Uuj 
and H.Friedrich, [Frlj . 

vi) In Section [4] the central part of our program is described and partially 
proved. In subsection 14. II an analogous result is proved for the Burger equation 
using, and somewhat extending, a previous result of S.Klainerman and one of the 
present authors (F.N.), [Kl-Nilj . Although the problem, in that case, is much 
simpler some of the basic ideas can be borrowed and transported. In the Burger 
equation case we prove that, due to the hyperbolicity of the equation, some 
apriori estimates hold for the Sobolev energy norms (with s — 2) up to a time T, 
depending only on these norms. Then it is proved that in the region of analiticity 
provided by the Cauchy-Kowalevski theorem it is possible to show that all the 
derivatives satisfy some appropriate estimates such that the series describing the 
analitic solutions have a convergence radius depending on the initial analitic data 

2 See for instance the detailed discussion in IPr-Rel . 

3 A slightly imprecise statement which, nevertheless, should be clear. 
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and on the Sobolev norms associated to the first derivatives, but independent 
from the point (of the analiticity region) around which we perform the series 
expansion. This is proved controlling, with a delicate inductive argument, the 
norms of all order derivatives. Once this result is achieved the analiticity region 
can be extended and as, again, the series convergence radius depends only on 
the initial data and on the Sobolev norms (bounded, via the apriori estimates, 
in terms of the Sobolev initial data norms) it can be proved that the procedure 
can be repeated in the whole time interval where the apriori estimates hold, 
obtaining the final result. 

In the second part of the section, subsection 14.21 we discuss how this approach 
can be implemented in the case of the Einstein equations. This requires a lot of 
technical work to which a subsequent paper is devoted. There we discuss in a 
very detailed way the problems one encounters and how we have solved them. 
We recall here some of these problems and give in subsection 14.21 an extended 
discussion of how they are faced. 

1) The first problem is the position of the initial data, in subsection 12.41 it was 
already discussed the constraints they have to satisfy, here we have to show 
how they can be given on the whole initial null hypersurface and not only in a 
portion of it, a generalization of what has been done in [Ca-Nilj . where we were 
interested only on Sobolev initial data. 

2) The local existence result has been proved in Section [3J the strategy to 
extend the analiticity region is to repeat the inductive mechanism which allows 
to control the norms of all derivatives again in an uniform way. This is the more 
complicated technical part; to achieve it we have to use the transport equations 
for the connection coefficients and the hyperbolicity of the Einstein equations, 
more precisely the a priori estimates for the integral norms of the Bel-Robinson 
tensor. The main lemma needed to prove our result, Lemma 14.71 is stated in 
subsection 14.21 while its proof and the subsequent steps to prove our result are 
written in the subsequent paper. 

To summarize this discussion we are convinced that to satisfy our goal of ex- 
tending as much as possible the analiticity region for the Einstein equations 
our gauge choice is the most convenient, even more as it seems so naturally 
intertwined to the characteristic problem. Moreover this formalism is perfectly 
suited to control the integral norms of the Bel-Robinson tensor and prove that, 
for small initial data, they can be bounded in the whole spacetime. 

Acknowledgments: The idea of trying to propagate as much as possible the 
analiticity of the solutions of the Einstein equations was suggested, years ago, 
to one of the author (F.N.) by S.Klainerman and was implemented together 
on the "toy model" of the Burger equation in [Kl-Nilj; moreover both authors 
are pleased to thank S.Klainerman one for the long and rich collaboration and 
(G.C.) for the useful period spent in the Math. Department of Princeton Uni- 
versity and for the long and profitable conversations he had with him. 



6 



2 The characteristic problem for the vacuum 
Einstein equations, assuming the spacetime 
foliated by outgoing and incoming null cones 

In this section and in the following one we present a way of writing the Einstein 
equations suited to study and solve the class of characteristic problems we are 
investigating. The basic idea is to assume that the spacetime we are constructing 
is foliated by a family of outgoing and incoming null cones, a foliation used in 
|K1-Ni2j which we believe very appropriate to study the characteristic problems 
for the Einstein equations. 

We will see that, with obvious differences, our approach is similar to the one 
used for the non characteristic problem when the spacetime is foliated by three 
dimensional hypersurfaces and the evolution part of the Einstein equations^ are 
first order equations in terms of the Riemannian metric of the hypersurfaces, 
gij , and their extrinsic curvature kij . 

Therefore the various steps required to accomplish our goal are in order 

i) Define the class of characteristic problems we are considering. 

ii) Define the gauge we use. 

iii) Identify the evolution equations in the coordinates associated to the 
gauge. 

iv) Identify the constraints equations in this formalism. 

v) the conservation of the constraints: Once steps i), ii), iii), iv) are clearly 
done, we show how the analytic solutions can be obtained in this characteristic 
case with the appropriate Cauchy-Kowalevski approach, that exactly as in non 
characteristic case we can define our analytic solution as a solution of a "reduced 
problem" and subsequently that, once the constraints are satisfied from the 
initial data, they are satisfied in all the existence region so that the analytic 
solution is really a solution of the (vacuum) Einstein equations. 

In the following subsections and in Section|3]we will concentrate on steps i),...v). 

2.1 The class of characteristic problems 

As we said in the introduction we are considering the case of initial data on 
a null hypersurface consisting of the union of a truncated outgoing null cone 
Co = C(Ao), see later for the "cone" definition and also |Ca-Nil] . and of a 
truncated incoming cone C = C_(vq) intersecting the previous one along an 
S 2 surface. Moreover we expect that analogous results can be easily obtained 
when the initial hypersurface is made by the intersection of two null hyperplanes. 
More delicate is when looking for the solutions of the Einstein equations inside 
an outgoing cone where we give the initial data. We believe that this problem, 
solved for the local existence in the previous work by by Y. Choquct-Bruhat, 
P.Y.Chrusciel, J.M. Martin-Garcia, quoted in the introduction, can be also faced 
with our technique, in our double null cone gauge. 

4 Therefore apart from the constraint equations. 
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2.2 The "double null foliation gauge" 



Let us recall what is the meaning of a "geometric" gauge choice in the non 
characteristic problem associated to a spacetime foliation. In that case, see also 
|K1-Ni2j , Lemma 1.3.2 , we can write the metric in the following way 

g(-, •) = -<P 2 dt 2 + g l3 (X l dt + dx^XHt + dx r ) (2.2) 

where <7y are the components of the riemannian metric induced on the generic 
spacelike hypersurface of the foliation. The unit normal to the generic hyper- 
surface is 

N= , 1 (l-X^) . (2.3) 

The coordinates {x 1 } ,i G {1,2,3} define a point p of a generic hypersurface 
labeled by a parameter r, p £ S T , and the diffeomorphism ^>n(S) associated to 
the vector field N send p to a point q G X r+ <5 with the same {a; 4 } coordinates. 
Up to now, although we have defined the coordinates {t, x 1 }, the "gauge" is not 
yet defined as we have not specified the lapse and the shift functions <J> and X. 
They are defined completely once that we define the hypersurfaces which foliate 
the spacetime. For instance in the choice used in |Ch-Klj . the hypersurfaces are 
assumed labeled by the time coordinate r = t and moreover they are assumed 
"maximal" which means that the trace of their second fundamental form is 
identically zero. The first statement implies that X = and the second one 
that the lapse function $ satisfies the hyperbolic equation 

A$ = \k\ 2 <f> . 

Similar considerations can be done for the "gauge" associated to the "CMC 
foliation" used by L.Andersson and V.Moncrief, |An-Mon"Tj . 

Following these ideas the "geometric" gauge we are choosing is associated to 
a double null cone foliation, see |K1-Ni2j and later on for its definition. The 
spacetime we are going to construct is foliated by a family of null outgoing 
cones C(A) and a family of null incoming cones C_{v) 1 more precisely a portion 
of null cones, see later on for details. We denote S(X, v) their intersections, 

S{\v) — C(X) nC(v) (2.4) 

which are two-dimensional surfaces diffeomorphicjfl to S 2 and with N and N_ 
the equivariant vector fields whose associated diffeomorphisms $ and $ send 
the S(X, v) surfaces to the analogous surfaces on the outgoing or the incoming 
cones respectively, 

${6)[S(X,v)\ = S(X,v + 5) 

§L(8)[S(\,v)] =S(X + 6,v) , (2.5) 

5 Here this is part of the assumptions we are doing about the foliation, later on we prove 
that this is true where the solution exists. 



8 



Let us also define the map ^(A, v) we will use extensively later on, 



(A, u):S 3p ^q = * (A, v)(p ) = *(u - n>)(£(A - A )(po) € S(A, i/) . (2.6) 

Once we have specified the foliation of the spacetime we can define two coordi- 
nates adapted to our gauge, namely the parameter A and v which determine the 
incoming or the outgoing cones. More precisely we denote these coordinates u 
and u and the outgoing and incoming cones C(A) and C_(y) are, respectively, 

C(X) = { P e M\u( P ) = X} , C(v) = { P eM\u( P ) = v} (2.7) 

where M denotes the spacetime. 

As the coordinates u, u allow to define the null cones it is important to know 
how to connect them to an arbitrary set of coordinates. This is provided by the 
eikonal equation, 

grdywduw = (2.8) 

where g^ v are the components of the inverse of the metric tensor written in 
an arbitrary set of coordinates x M . Then u = u(x) and u = u{x) are solution 
of the eikonal equations with appropriate initial datc[^ such that their level 
hypersurfaces are, at least locally, the null cones of our foliation. It is well 
known, see for instance |K1-Ni2j . that the vector fields 

are the tangent vector field of the null geodesies generating the outgoing and 
the incoming cones, satisfying 

D L L = , D k L = . 

Their scalar product defines the scalar function fi, 

g(L,L) = -(2n)- 2 (2.9) 

and the vector fields e$, e^, 

e 3 = 2QL , e 4 = 2QL , (2.10) 

satisfy the relation g(e3, e 4 ) = —2. Finally the equivariant vector fields previ- 
ously introduced, N, Njj have the following expressions 

N = fie 4 , N = Qe 3 . (2.11) 

Once introduced the coordinates u,u we are left with defining the remaining two 
coordinates, we have to interpret as angular coordinates, which allow to specify 



6 Different choice of the "initial data" for l2.8l on the external outgoing and incoming cones 
give rise to different null cone foliations. 

7 It is immediately to check that they satisfy eqs. 12.51 
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a point on each surface S(X, v). The procedure is similar to the one envisaged 
previously in the non characteristic case, to define them we have to define a 
map which sends a point p of "angular coordinates" {uj a } , a G {1, 2} belonging 
to the intersection of the "initial cones" S(Xq,vo) = C(Ao) flC(i/o) to a point 
g £ S(X, v) to which the same coordinates are assigned. The map we choose is 
made in the following way: first we move, starting from a point p £ S(Xo, vq) 
of angular coordinates {oj a }, along the C_{vq) cone using the integral curves of 
the vector field JV up to a point q' £ S(X, vq), then applying the diffeomorphism 
generated by the vector field N we move "inside" the region of the spacetime up 
to the point q £ (A, v) and assigne to this point the same angular coordinates 
{w a }; formally, for any a, 

w°(*(A, v){p)) = uj a {<$>{v - ^ )$(A - Ao)(p)) = w a (p) ■ (2.12) 

Once we have introduced the coordinates relative to our gauge we have still to 
write, in these coordinates, an explicit expression for the metric tensor. This 
can be done defining a null moving frame, {e4,e3,eyt}, A £ {1,2} adapted to 
this double null foliation where 

e A = e a A -^, A€{1,2} , (2.13) 

are orthonormal vector fields tangent at each point p £ M. to the surface S(X, v) 
containing p and , are null vector fields orthogonal to the e^'s, outgoing 
and incoming, respectively, which, basically, means "tangent" to the null hyper- 
surfaces C(A) and C(^)o Moreover, recalling the meaning of the vector fields 
N and N_, it follows that 

The explicit expression of N_, and therefore of is somewhat different as in a 
curved spacetime these two vector fields do not commute; it can be proved, see 
|K1-Ni2j that N_ must have the following expression 

* =i + xa i and therefore 63 = h {I + xa i^) • (2 - 15) 

where the property and the equation that the vector field X has to satisfy will 
be discussed later onf^ 



Once we have the explicit expression of the null frame in the adapted 10 ! coordi- 
nates we can write the metric tensor in these coordinates obtaining 

g = ~2fl 2 (dudu + dudu)+-/ ab (X a du + duj a )(X b du + duj b ) (2.16) 



8 e4 is, at same time, tangent and normal to C(A) and analogously ez with respect to Cjy). 
9 The definition of the coordinates ui a , eq. 12.121 implies that X = on C_ . 
10 Here with adapted we mean both adapted to the "gauge" we are defining, both to the 
leaves S(A, v) of the null cones. 
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where j a b are the components of the induced metric on S(X, v). 

Remark: Observe that this is not the more general metric we can write "adapted 
to the foliation, in fact there are only seven metric components different from 
zero. This follows as we have chosen the coordinate u such that the N vector 
field was . We could nevertheless define the coordinates in such a way that 

N — + Y a , in this case repeating the previous argument the metric turns 
out to be, written in the same u,u coordinates, 

g = -2n 2 (dudu+dudu)+ lab {X a du + Y a du + duj a )(X b du + Y b du + duj b ) 

with ten components different from zero^^ This last expression is the one anal- 
ogous to the metric [2J% written for the non characteristic problem 

g(., .) = -$ 2 di 2 + gijiX^t + dx l )(X^dt + dx j ) 

while the previous expression \2.16\ corresponds to the one used for instance in 
]Ch-KHI associated to the "maximal" foliation, namely 

g(., .) = -$ 2 dt 2 + gij dx l dx 3 . 

Comparing the metric expression 12.161 with the one in 12.21 we see that il plays 
basically the role of the lapse function and X of the shift vector. Exactly as 
in that case we expect that the gauge choice is completed once we are able to 
determine these functions. In the non characteristic case this could be done in 
different ways, for instance in |K1-Ni2j the shift vector was imposed equal to zero 
and the lapse function had to satisfy an elliptic equation. In other cases, see for 
instance jCh-Klj or jAn-Monlj l 12 ! both quantities had to satisfy some differential 
equations. In the present characteristic case we will see that both £7 and X a 
have to satisfy some first order differential equations. We will obtain them in 
the next sections after we write the Einstein equations in these coordinates. 
To do it we use in a systematic way the structure equations for a Lorentzian 
manifold. 

2.3 The Einstein equations in the {A, v, Q, cp} coordinates 

In this section and in the following section we fulfill the first goal of this paper, 
namely we write the Einstein equations in a way suited to the class of character- 
istic problems we are considering and connected to the gauge choice discussed in 
the previous subsection. We recall shortly the main properties of our approach, 
i) As we said before we assume the spacetime foliated by null cones, outgo- 
ing and incoming as prescribed by the chosen gauge. Therefore the Einstein 
equations become evolution equations for the metric components j a t of the two 

11 Nevertheless not all independent due to the nature of the double null cone foliation. 

12 In their case the spacetime is spatially compact and the reduced equations are the evo- 
lution equations (2.5a), (2.5b) together with the elliptic equations (2.8a), (2.8b) relative to 
N and X which specify the gauge; given a solution of this set, the constraint equations are 
proved to hold everywhere in the spacetime once they are assumed for the initial data. 
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dimensional leaves S(X, v) which foliate the null cones. Moreover as in the non 
characteristic case the choice of the gauge will imply that the lapse function, O, 
and the shift vector X have to satisfy some equations. 

ii) As in the analogous formulation of the non characteristic case we express the 
Einstein equations as a system of first order equations; we define the family of 
first order equations which express the Einstein equations and separate them in 
two groups, in the first one we collect those equations which can be interpreted 
as the evolution part of the Einstein equation^] and in the second one those 
which can be interpreted as constraint equations. This is the central goal of this 
subsection. 

iii) The third aspect to remark is that we write our equations in terms of the 
coordinates u = A, u = v, cj a , where v and A are, as defined in the previ- 
ous subsection, affine parameters for the null geodesies along the null outgoi ng 
and incoming cones and we do not use the more standard time coordinate, t ri 
With this choice our equations have a more geometric "flavour" as all the quan- 
tities we introduce are connected in a direct way to the geometric properties 
of our foliation and the equations are transport equations along the null cones. 
This approach and this formalism turn out to be very appropriate, as already 
discussed in |K1-Ni2j . to obtain apriori estimates for "energy- type" quantities. 
This, moreover, will allow us to prove the main goal of this work, namely, 
as discussed in the introduction, that the analytic solution of the characteris- 
tic problem has a much larger existence region depending on the "hyperbolic" 
apriori estimates, a region which can be unbounded if the initial data have H s 
norms, with appropriate s, sufficiently small. 

To write the Einstein equations in a way satisfying i) ii) and iii) we use in- 
tensively the structure equations, see for instance, |Sp| , Vol 2, adapted to a 
Lorentzian manifold. 

2.3.1 The structure equations 

We recall some general aspects of the structure equations. We denote the null 
orthonormal frame in the following way: 

{e( Q )} = {e a } = {e (1) , e (2) , e (3) , e (4) } (2.17) 

where, 

e (3 ) = 2QL, e (4 ) = 2QL , e (1 ) = e d , e (2 ) = ( 2 -18) 

and 



= {o<*} = {9W,ew,ew,oW} (2.19) 



13 Together with the equations for and X. 

14 The simplest analogy is solving the two dimensional homogeneous linear wave equation 
written as: 

d 2 

u = . 

d\v 
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are the corresponding forms satisfying 

6< a He w )=8% (2.20) 

and it follows that 

= -\g^ v , ol 4) = -\g^ v . (2.2i) 



We define 



0^7 

R(e Q e^)e 7 = R 7Q(3 e 5 (2.22) 

where D is the connection of the spacetime associated to the Lorentz metric g, 
D eQ e^ is the covariant derivative of the vector field ep in the direction e a and 
R is the Riemann tensor (here the first greek letters are "names" and do not 
denote components) 

R(e ce e /3 )e 7 = D ea (D e(3 e 7 ) - D e|3 (D eQ e 7 ) - D[ e ^ e)3 ]e 7 

Defining the following one and two forms: 

"1 = V 

n% = iR^Afl 4 (2.23) 

we have the following result, see [Sp], Vol.2, whose proof is in the appendix to 
this section. 

Proposition 2.1 cj^ and £lp satisfy the following structure equations 

d6 a = A (9 7 (2.24) 
duj* = -u' a A u° + 0* , (2.25) 

called, respectively, the first and the second structure equations. 

The knowledge of a null orthonormal frame in a whole region is equivalent to 
knowing the metric in that region, therefore the first set of structure equations 
can be thought as "first order equations" for the metric components. Viceversa 
the one forms ujp are connected to the first derivatives of the (components of 
the) moving frame and, therefore play the role of the first derivatives of the 
metric components; the second set of structure equations represent first order 
equations for these first derivatives F^l 



15 Obviously choosing a coordinate basis {g^rr, g~T> Jjfsi eS 3 '^ ^ e defined in 12.221 are 
just the usual Christoffel symbols, T„ p , D a = D p is the covariant derivative with respect 

to d„ and R* Q/3 = KV e 7 e £ e F£ ' 

16 Observe that the second group of structure equations depends, through f22, also on the 
Riemann tensor components. This could suggest that expressing these equations as partial 
differential equations for the various components, to have a closed system of equations one 
should also consider the Bianchi equations for the Riemann tensor and the Riemann compo- 
nents as independent variables. Although this could be done, see for instance |Frl| . this is not 
what we do, as we discuss in great detail later on. 
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To have an explicit expression for the structure equations we recall the definition 
of the "connection coefficients" (sometimes called "Ricci coefficients"). In the 
defined null orthonormal frame they have the following expressions: 

Xab = g(D eA e 4 ,e B ) , x ab = g(D eA e 3 ,e B ) 

i A = ^g(D e3 e 3 ,eA) , U = ig(D e4 e 4 ,e A ) 



w = ^g(D e3 e 3,e 4 ) , w = ig(D e4 e 4 ,e 3 ) 
1 , > 1 



(a = ^g(D eA e 4 ,e 3 ) 

They are 2-covariant tensors, vectors and scalar functions defined on the two 
dimensional surfaces S(X, v). 

In terms of these quantities the one forms, u^, have the following expressions: 

wf = , u>i = -2^ + 2cj8 3 - ( a 9 A , < = V A 8 4 + lx BA B 

u 3 = 2lo6 4 - 2lj9 3 + ( a 6 A , u 3 = , ui\ = VA 9 3 + ^ X baO b 

lu a = 2 Va 9 3 + X baO b , <O3=2v A 0i + X B A 0B ( 2 ' 26 ) 

The one forms oj^ have a different expression which do not depend on the 
connection coefficients we have introduced, 

lo a = g(p 4 e B , e A )B 4 + g{p 3 e B} e A )0 3 + g(D c e B , e A )0 c . (2.27) 

It is a long, but simple task to write the structure equations in terms of the 
metric components and the connection coefficients. The first set of structure 
equations 



becomes: 



d6 a (e p ,e 7 ) = -co% A6 s (ep,e 7 ) 



{a,/3, 7 }-{A,S,4}: — lab = 2ft Xab 

{a, /3, 7} = { A, B, 3} : — 7ab + C xlab = 2% ab 

{a, (3, 7} = {A 4, 3}: |-X C = 4fl 2 Cce^ (2.28) 

{a, /3, 7} = {4, 4, 3}: logfi = -2Qu, 

{a, /?, 7} = {3, 3, 4} : f A + X a ^A log = -2ftu; 
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{a, A 7} 
{a, ft 7} 
{a, A 7} 
{a, ft 7} 



{4, A 4} 
{3, A3} 
{4, A3} 
{3, A 4} 




a = . 



(2.29) 



Remark: To obtain equations \ 2.28\ and \ 2.29\ we have used only the relation 
63 = 2f2L , 64 = 2QL. Recall that the structure equations for a moving (null 
orthonormal) frame do not imply that {e^, e^}, A S {1, 2} or {e^, e^}, A € {1, 2} 
are integrable distributions. In our present case, with the previous definition of 
e3, e4, this is true. 

As our purpose is to write the Einstein equations as a system of first order equa- 
tions for the component of the metric and their derivatives, equations l2.28l of the 
first set of structure equations tell us exactly that we have to look for first order 
partial differential equations for the connection coefficients {x, X, C, w, w}. They 
correspond in fact to the derivatives along the inward (or outward) direction 
with respect to the null cones outgoing or incoming, of the metric components 
7a6, of the lapse function and of the shift vector. Once we have written the ap- 
propriate partial differential equations for these connection coefficients we will 
show that we reduce to a closed set of equations in the metric components and 
in these connection coefficients. To find the appropriate equations we rely on 
the second set of the structure equations, 



These equations can be written in a more explicit way in terms of the connection 
coefficients. One has, nevertheless, to remember that the structure equations 
are identities, valid in a generic manifold with a (Riemannian or Lorentzian) 
metric. They can also be seen as "integrability conditions" for the existence of 
a moving frame in the whole manifold. One has to remark that in their right 
hand side the terms f2*(e a ,e^) appear which can be written in terms of the 
various components of the Riemann tensor. Therefore these equations can be 
seen as identities defining the Riemann tensor in terms of first derivatives of the 
connection coefficients, but if we impose some condition on the Riemann tensor 
they become partial differential equations to solve with respect to the connection 
coefficients. The condition to impose to the Riemann tensor, or more precisely 
to f2*, are that the vacuum Einstein equations have to be satisfied, namely that 
the corresponding Ricci tensor is identically zero. Therefore we have to look to 
the explicit expression of these equations under this condition. 
It is a long, but standard and certainly not ne to realize that the structure 
equations have the expressions written in the following where we indicate also 

17 What is certainly more uncommon is that we do not use a subset of these equations to 
get some norm estimates for various terms assuming we already have a solution, but to solve 
them, which requires a delicate choice of the subset. 



(du4 + w% A u°)(e a , ep) = Q S Je a ,e ) . 
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the Ricci components (= 0) to which they are associated^ 



[R(e 4 ,e 4 )=0 :] 

D 4 tr X + i(trx) 2 + 2wtr X + \x? = 



[R(e 4 ,e A ) = :] 
p 4 ( + Cx + trxC 



divx + ytrx + PJ7 log Q = 



[S A B'R-(eA,eB) = 



D 4 trx + -trxtrx 



2wtrx + x-X + 2div(C-yiog^)-2|C-yiog^| 2 = 2p 



[R(e A ,e B ) = :] 

P4X+ \^XX+ ^trxx - 2^x+ f%(C~ f^og 0) - (C - ^log fi)§(C~ ?log O) = 



[R(e 3 ,e 4 )=0 :] 

D 4 w-2^-C ■ yiogO-||C| 2 + ||yiogO| 2 + I(K + itrxtrx-^X • X) =0 
[R(e 3 ,e 3 )=0:] 



[R(e 3 ,e j4 ) = 0:] 

PaC + trxC + CX + 4ivx - ^tr* - P 3 ? log O = 
[^ B R(eA,es) = :] 

p 3 trx+trxtrx-2a;trx-2o}iv(C + yiogO)-2|C + yiogO| 2 + 2K = 



where D3 = D e3 ,D 4 = D e4 , JZ)3,JZ) 4 their projection on the tangent spaces 
TS(X,v), J/ the covariant derivatives associated to the metric 7 induced by 

18 Remind that we used the relations between r), r\ and £ , the expression of oj and in terms 
of Q and the fact that § = £ = 0, obtained from the first set of structure equations, see for 
instance Kl-Ni2 , Chapter 3. 



D 3 trx + ^(^X) 2 + 2wtrx + |x| 2 







(2.30) 
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g on the surfaces S(X, v), K is the curvature of these S(X, v) surfaces and 

R(eA, e B ) = R(eA, e B ) - 2~ 1 5 A B'R{eA, e B ) ■ 

In the appendix we write the general form of the second set of structure equa- 
tions and show how from them, the first set and our gauge choice equations [2301 
follow. 

We have now the explicit expression for the first set and the second set of the 
structure equations, namely 12.281 and 12.301 They are 24 equations, ten from the 
first set and fourteen from the second one while we have only sixteen unknown 
functions f2, X, 7, x, X? d u i The fact that there are more equations that un- 
known functions is not, in the present real difficulty as the structure 
equations are automatically satisfied in a Lorentzian manifold and, therefore, 
also in a vacuum Einstein manifold. Therefore we have to choose a subset be- 
tween them which forms a complete set of equations for the sixteen unknown 
functions and then prove, as expected, that the remaining equations play the 
analogous role of the standard constraint equations and are automatically sat- 
isfied once they are imposed on the initial data. The nature of the constraint 
equations is discussed in detail in subsection 12.41 The sixteen equations we 
choose are: 

|^-2ftx + £x7 = 
dl °d\ l + lo & n + 2n ^ = 

p 3 trx+trxtrx-2wtrx-2chV(C + yiogfi)-2|C + 'y 7 logO| 2 + 2K = 
^3X + ^trxx+itrxx-2wx-y®(C+?log^)-(C+^ogri)§(C+yiogr!) = 
PsC + trxC + Cx + divx - ytrx - P 3 fiog Q = (2.31) 
D 3W -2w W + C ■ yiogO-||C| 2 + ||yiogO| 2 + I(K + itrxtr X -^X • *)=0 



— + m 2 z = 

ov 

^) 4 trx + trxtrx - 2wtrx + 24iv(C-yiogfi)-2|C-yiogfi| 2 +2K = (2.32) 
P4X + 5 tr XX + - trxx - 2^X + T® (C - ^log 0) - (C - flog SI) § (C - ^log fi) = 
D 4 ^-2c^-C • yiog^-^|C| 2 + i|yiogfi| 2 + i(K + itrxtrX"^x • X) =0 . 
Remarks: 

i) The equations 1 2. 3 1\ and \2.32\ are appropriate, as we will see in the following, 
to apply Cauchy-Kowalevski theorem and find analytic solutions. To write them 
as first order p.d.e. equations for the tensor components requires still some more 
work due to the presence of the Gauss curvature K which depends on the second 
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angular variables of the metric 7. This will be discussed in the next subsection. 



ii) The second important remark is that once we have solved equations ] 2. 3 1\ and 
\2.32l looking at \2.30\ our analytic solution is such that 

R(eA, e B ) = , R(e j4 , e 4 ) = , R(e 3 , e 4 ) = 

and we have still to prove that the remaining Ricci equations are satisfied. This is 
discussed in detail in subsection \2.4\ where we show that the remaining equations 
to be satisfied have to be considered as "constraint equations". 

Hi) Observe that in this approach there are ten independent connection coeffi- 
cients, x, They, basically, correspond to the second fundamental form 
kij of the "maximal foliation gauge" or of the "CMC foliation gauge". The dif- 
ference is that k has only six components. This is due to the fact that in those 
cases the foliation is made by only one family of hypersurfaces St while here 
there are both the C{\) and the C_{y) null hypersurfaces. If we consider only 
the {C(A)} foliation, x, are the S-tensors corresponding to k (x, (,u in the 
opposite case). In both situations there are six components, as expected. 

Equations l2.31[[2.32l are perfectly defined as tensorial equations, but, to consider 
them as p.d.e. equations, they have to be written as equations for the tensor 
components. In this case they do not maintain exactly the same expressions. 
There are many ways to rewrite these equations as standard partial differential 
equations whose unknown are the components of the various tensors involved, 
for instance one could choose a Fermi transported null orthonormal framj^l as 
was done in |K1-Ni2] ; here we present a more general approach using the diffeo- 
morphism ^(A, v) previoulsly introduced to map, via the pullback associated 
to ^(A, v), these equations on a manifold So x R 2 where the equations become 
equations for the various components in the angular variables and in the vari- 
ables A and v, which are just the parameters of the diffcomprphism ^(A, v). 
Here we state the result and its detailed proof is in the appendix. 

Proposition 2.2 In the coordinates {A, v, 6, <f\ associated to the "double null 
foliation gauge", eauations \2.31\ and \2.32\ written for the various of metric and 
connection coefficients components have the following expression: 

+ fx log^ + 20 lu = 
— -p + Otrxtrx - 20wtrx + f x trx ~ 20o}rv(C + y io S n ) ~ 2 ^IC + yi°gO| 2 + 20K 

|| - + + ^|^x + (f x log fi)* - n(x • xh - n f®(( + f^g n) 
-r>(c + yio g n)g(c + ^iogO) + £ x x=o (2. 



19 Nevertheless this is possible only with respect to a null direction, but not simultaneously 
to both. 
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dC dWlos ft 

^ + trxC + Vtyvx - nftrx - ^ + Qf log fl-x_ + CxC - CxflogVt = 

^+f x u-2nuLo~n\(\ 2 +ncyiogn+^n\fio g n\ 2 +^n(K+^^^ 



^ + m 2 z = o 

Of 

dtrv „ 

+ Otrxtrx - 2ftwtrv_ - 2f2d4v(-C + f log ft) - 2fi| - C + flogVt\ 2 + 2OK = 

dx iltrv fitrx dlogfi 

£ -^x + -^x + -J-x - m ■ xb ny®(-c + f log fi) 

-fi(-C + yiogO)®(-C + yiog^) = (2.34) 
^= - 2ttc^- -tt|C| 2 - fX • yiogtt + -ft|yiogtt| 2 + -ftf K + -trxtr X --* • x] = . 

Equations l2. 331 [2.341 are not yet a system of first order equations. In fact K, the 
curvature of the two dimensional surfaces S(X, v), depends on second tangential 
derivatives of 7. Moreover in the (transport) evolution equations along C_{y) for 
t r X; Xi C an d in the (transport) evolution equations along C(A) for tr% and x 
the second derivatives of log f2 with respect to the angular variables are present. 
To have a real first order system we define some new independent variables and 
their evolution equations, namely: 

«... = 0/y(v) , w..=p.X(-) , ^. = ^.logfi • (2.35) 

Their evolution equations are obtained deriving the evolution equations of / y(>, •), 
X(-) and logfi, the second one in the outgoing direction, the other two in the 
incoming one. The unknown function w... is introduced as in the evolution 
equation for v the second tangential derivatives of X appear. These equations 
do not contain more than first derivatives of the previous unknown variables 
and this transforms the system of equations into a larger system of first order 
equations. It is a matter of computation, which we report in the appendix to 
this section, to obtain the following evolution equations for ip ai v ac ib and for w a b'. 

^-Va = ~2Qd a W - 2Qu^ a - {f a Xf^ c - X C f c lJj a 



-Q^ v cba = -(d c X d )v dab - d x v C ab + d c w ab + d c w ba + 2fld c x ab + 20-0^2^36) 

= -8n 2 v«a - m 2 d a ( b + 2n^ aXbc x c + 2n(d aXb c)x c . 

ov 

We write now the final system of first order equations for the various tensors 
components omitting the indices to simplify the notations, 

d 1 diogn ox 

v = , — ip = Q , — w = 

aw ouj ou! 
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d log VL , . 

B +ip(x) + 2nu = o 

OA 

dv 

d\ 
dip 
d\ 
dtvx 
~dY 



+ fx v + $X)-v- S($ ®w) - 2fl@ ®x - 2fl ip(Z>x = (2. 
+ f x ^ + 2Qun/j + tp{fX) + 2Q,$ui = 

Qtrxtrx - 2ftwtrx + f x ^X ~ 2fid4v(C + ip) - 2Q\( + ijj\ 2 + 2ftK = 



d\ ritrx fitry 

-± + c xX - -^= X + ~y^x - u x - fi(* • xh - n ?®(C + tf) - fi(C + V0®(C + ^) 

<9C 1 

+ + fitrxC + ^4^vx - Tr^trx +2Qunp + 2f^w + fi?/"X = 

^ + ^ W -2^ W -^|C|V^C> + ^l^| 2 + ^(K + itrxtrx-ix-x) = 
av 

+ SO^C + 4fi 2 ^®C ~ 20^(g)(x-X) - 2n($®x)-x = 
9trx n 

-^-= + fitrxtrx - 2fljtrx + 2ftqMvC-2rJc}ivV> - 2fi|C-^l + 2ftK = (2. 
d\ fitrv f2try ~ ~ 

# - -A + -t^x - 2ftwx - m ■ xh + ny<8(c-^ - a C-V> ® C-VO =o 

01/ 2 — 2 — — 

where X is the covariant vector X a = j a bX b , S means symmetrization, V®W 
is twice the traceless part of the symmetric tensorial product S(V <E> W), K has 
to be thought as a function of 7, v and SkiW\ 



2.4 The first order system of equations as solutions of 
the vacuum Einstein characteristic problem, the con- 
straint problem. 

The first order system of p.d.e. equations 12.371 12.381 describes a characteristic 
problem which can be solved via the Cauchy-Kowalevski theorem (its charac- 
teristic version, as discussed later on in Section [3]), giving the inital data on the 
two null hypersurfaces Co and CpFI As we said before, see remark ii) after 
equations I2.32[ the equations we want to solve are not all the equations asso- 

20 <j) is the ordinary partial derivative with respect to the angular variables, u) a , and y is 
the Levi-Civita connection with respect to 7. 

21 Which initial data can be given freely and which constrained is a delicate point we discuss 
in detail in subsection 14.21 
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ciated to R M „ = 0. Therefore we have to determine under which conditions a 
solution of the equations 12.371 l2~38l is a solution of the Einstein equations. Ob- 
serve that, looking at the structure equations, apart from equations 12.371 [2.381 
the following equations have to be satisfied by the vacuum Einstein equations: 

£-2O* = , ^ + 2^ = 
Ov ov 

dC 

R(e 4 ,e A ) = : + QtrxC - ^d/ivy + OTtrx - 2nuJip - 2Vlfuj - Qip-x = 
ov 

R(e 4 ,e 4 )=0: ^ + ^tr X + 2^tr X + «|x| 2 = (2.39) 
Ov 2 

<9trv i7try 

R(e 3 ,e 3 )=0 : ^ + ^try + ^tr* + 2fWx + ft|x| 2 =0 . 

The first two equations of I2.39I are at the "level" of equations for the metric 
components and follow from the first set of the structure equations. This is 
somewhat analogous to what happens in the maximal foliation gauge used by 
D.Christodoulou and S.Klainerman, |Ch-Kl) . where, once we impose trfc = on 
So, one has to prove that trfc remains equal zero on any t-constant hypersurface, 
justifying the definition of the maximal foliation gauge. In other words proving 
that the first two equations, once satisfied on Co, are satisfied on any C(A) 
shows that we are in the double null foliation gauge. 

The remaining three equations have to be satisfied to make the components of 
the Ricci tensor, R(e4,e y i), R(e 4 ,e 4 ) and R(e3,e3), identically zero. Here it 
is appropriate to introduce the notion of "signature" for the various functions 
involvedri 

Definition 2.1 We call "signature" of the various connection coefficients the 
number of times the null vector e 4 appears in their definition minus the number 
of times e 3 is present. Each derivative along Cq increases the signature by one 
and viceversa for each derivative along C . 



Observe that the last three equations in 12.391 are at the level of connection co- 
efficients and have signature +1, +2 and —2. This can be interpreted as the 
indication that in these equations there are no "derivatives" with respect to 
the transverse directions, e 3 for C(A) and e 4 for Civ). These equations have, 
therefore, to be seen as constraint equations and we have to prove that, if sat- 
isfied from the initial data on Co and C , they are satisfied on each outgoing or 
incoming cone, respectively^! This is the content of the following lemma which 
connects the solutions of 12.371 T2.38I to the solutions of the Einstein equations, 

Lemma 2.1 Let * = (j ab , log O, X a , v c , abl w ba , \ab, Ca, w,X af) ,w) be a solution 
of the first order system made by equations \2.37\ \2.38[ 

22 This was introduced by D.Christodoulou and S.Klainerman, Ch-Kj] . for the null compo- 
nents of the Riemann tensor. 

23 Observe that the effect of a coordinate choice and of the choice of the system of equations 



make the set of equations [Z39] asymmetric with respect to the A, v interchange. It is also easy 
to see that we have a certain arbitrariness in choosing the first order system, for instance one 
could interchange the role of the 'V and "A" directions. 
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If {"fab, log trx, £} are a solution of the first four equations of \2.39\ on Co and 
tr\ is a solution of the last one on C_ , it follows that they are solutions of the 
same equations on any cone C(A) and C[y) respectively. 

Proof: We show that there exist first order transport equations along A for the 
left hand sides of the first four equations of 12.391 and a transport equation along 
v for the left hand side of the last equation. Therefore, if these expressions are 
zero on Co or on C , they are identically zero for all A, v values. A way to obtain 
these transport equations is just a long computation using equations 12.371 12.381 
We write only the proof for the second equation of the first line of 12.391 

d /Slog n nn \ d fdlo g n\ /<91ogft\ Buj , n ln . 

ox {-w- + m i = *, Hr ) + 2 Hr r +m ox (2 - 40) 

Using again equations 12.371 12.381 the terms in the right hand side become 
d ( <91ogf2\ d , „^ _ , _ x 

1 — (-2fiw-a x iogfi) 



dv V dX J dv 



± log n W 2n^ - d -fdj og n-dj^l) (2.4i) 

ov J ov ov \ ov J 

2 ( dl Q g ^ ^j ^ = 2ttuj{-2nu - d x \og fi) = -Atfuju - 2fiw9 x log Q, . 

Substituting in [2222 and denoting I = (^§f^ + 2£Iuj\ we obtain 



dX 

— = -Wyl - 2Q,oj1 + 2n 
dX rx - 



(2.42) 



and, using again equations 12.371 12.381 the term in the square bracket is identi- 
cally zero. Therefore I satisfies the following equation 



ar 

dX 



y x l + 2VLul=Q (2.43) 



which implies that I is equal to zero on every outgoing cone C(A) provided it 
is set equal zero on Cq. 

To complete the proof of the lemma at the level of the metric components, let 
us consider the first equation of 12.391 Proceeding as before and denoting 

t- dlab on 

-Lab — ^"Xab , 

OV 

a long, but straightforward computation shows that I a }, satisfies the equation: 

— ^ - QtrxZab + Cxlab = (2.44) 

which implies again that if ^Sjff- — 2f2% & = on Co then this relation holds on 
any C(A). To complete Lemma [2.11 the same result has to be proved for the 
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remaining equations in 12.391 If we proceed as before the computation will turn 
out very long and laborious. This can be avoided observing that this result 
follows by a straitghforward application of the Bianchi equations. In fact let us 
consider the Lorentzian manifold with metric 

g(-, ■) = \X\ 2 d\ 2 - 2Q 2 (dXdu + dvdX) - X a (d\duj a + dcu a d\) + ~/ ab dcu a du b , 

where X, fi, 7 satisfies equations [533 [23H From these equations and the result 
just stated it follows that the components of {Xab, Ca, w, x ab , w}, can be 
interpreted as the connection coefficients associated to this metric. Therefore, 
as they satisfy equations 12.371 12.381 it follows that, see equations 12. 30i the null 
components R(e j 4, eg), R(e3, e^) and R(e3, &a) of the Ricci tensor are identically 
zero. To prove the remaining part of Lemma 12.11 amounts to prove that also 
R(e4,e4), R(e J 4,e4) and R(e3,e3) are identically zero, provided they are equal 
to zero on the initial hypersurface. To prove this result we use the contracted 
Bianchi equations. In fact from them one deduces the following identities 

D»R^ - \d v R = . (2.45) 

Denoting {e3, e^, ca}, A 6 {1, 2}, a null orthonormal frame and writing 

A 

equation 12.451 can be written, mutiplying it with and es respectively, 
-l(D 3 R^)e%eZ + Y,(D A R^)e^i ~ l(D*R^)e*<& = 

A 

-\{pjt,J)4^ B - \{D z R^)e^ B + J2(D A R^)e1eB + \{D sR^e^ 

A 

-Y.V D b R ^a^a=^ ■ (2-46) 

A L 

Rewriting these equations as transport equations for the various null Ricci com- 
ponents, from the first set of structure equations for the null frame, see for 
instance [Kl-Ni2] . Chapter 3: 





1 

= J A e B + ^XABe 3 


1 


Dac 3 


= X AB e B + de 3 , 


T>Ae± = XABeB - Ca&a 


D 3 eA 


= p^e A + r] A e 3 , 


D 4 e^ = p4e A +V A e i 


D3e 3 


= (D 3 logO)e 3 , 


D 3 e 4 = -(D 3 logft)e 4 + 2r/ B e B 


D 4 e4 


= (D 4 logrj)e 4 , 


D4e 3 = -(D 4 log^)e 3 + 2r/ B e B 



and recalling that all the null Ricci components with signature —1 and are 
already equal to zero, equations 12.461 become: 

- + <9xJR(e 4 ,e 4 ) - trxR(e 4 ,e 4 ) + 4wR(e 4 ,e 4 ) + 2f A R(e Al e i ) 
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2(f?s - g(f A e A ,e B ) + y s logQ)R(e B ,e 4 ) = (2.47) 
d 



. jX . 9 x JR(e s ,e 4 ) - trxR(eA,e 4 ) - x(eA,e B )R(es,e 4 ) - 4wR(e B ,e 4 ) 

+5(^36.4, e B )R(e_B, e 4 ) =0 . 

From the second equation and the assumed initial conditions it follows that 
R(e B , e 4 ) = 0, which, substituted in the first equation, implies that R(e 4 , e 4 ) = 
0. The proof that R(e 3 , 63) = goes in the same way, is somewhat simpler and 
we do not report here. Once Lemma |2"7T1 is proved it follows that $ is a solution 
of the vacuum Einstein equations. 

Remark: The previous discussion makes cristal clear, in the Einstein equations 
characteristic problem, which are the equations we have to consider as evolution 
equations and which have to be interpreted as constraint equations, which is 
enough to satisfy on the "initial data" to have them satisfied everywhere, The 
first ones are equations \2. 31\\2. 38\ while the "constraint equations" are equations 
\2.39\ which do not involve inward (outward) derivatives of the initial data. 
One has also to remark that if we consider only the equations \2.37\ \2.38\ and 
we do not care about initial data satisfvina \2.39\, we are still considering a well 
defined characteristic problem whose solutions nevertheless do not define an Ein- 
stein vacuum spacetime. Nevertheless as this is a characteristic problem in itself, 
even in this more general case the initial data have to satisfy some constraints, 
namely the initial data associated to {7, ft, v, ip, X, C> w } have to satisfy the con- 
straints prescribed by eqs \ 2.31\ on C = C_(vq), while they are given in a free 
way on Co = C(\q) and the opposite has to be imposed for the initial data of 
{X,w,x,ul}. 

In conclusion one has to recognize that, in some sense, the characteristic problem 
for the Einstein equations has two kind of constraint equations that the initial 
data have to satisfy, the first one connected to the more general problem 12.371 
12.381 and the second one to the requirement that also equations 12.391 have to be 
satisfied. We summarize this discussion in the following theorem, 

Theorem 2.1 Let * = (j ab , log ft, X a , v Ctab , ip a , w ba , Xab, (a, w, X ab i w) be a so- 
lution, in a region, {(\,v)\{\,v) G [0,A] x [0,17]}, of the characteristic first order 
Cauchy problem made by equations \2.37\ \2.38\ with the initial data on the null 
hypersurface S = Co U C_ satisfying on S, beside equations \2.37\ \2.38[ con- 
sidered as equations on C_ and Cq respectively, the constraint equations, see 

d'y dlos Q 

On C : ^ - 2Q X = , + 2fiw = 

ov ov 

dtrx , ffirx,„, , na-i _, , ok -,2. 



dv 2 



^ + ntrx(-nq v ivx + n@tr X + d @ l ° gn -nffio g o-x=o 

ov ov 
dtrv fltrv 

OnC : -^ + ^trx + d x trx + 2n^trx + n\x\ 2 =0 , (2.48) 
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then in the same region the metric tensor 

g (., •) = \X\ 2 d\ 2 - 2Q 2 (d\dv + dvdX) - X a (d\doj a + duj a d\) + j ab duj a duj b (2.49) 
is a solution of the Einstein vacuum equations. 

Summarizing we collect here all the constraint equations the initial data have 
to satisfy, namely equations I2.37| 12.381 and 12.391 

On C : 

$ 1 -v = o , fix -w = o , ^iog n - ip = o 

£-2nx = , ^+4O 2 Z = , ^ + 2^ = 
ov ov ov 

dv 2 



-trx + 2r2wtrx + ^|x| =0 



+ OtrxC ~ ^<Mvx + Ofex + ^ log ° - n^logn.v = (2.50) 
w ' Ov 



<9z/ 



fitrxtrx - 2Slo;trx + 212o}lvC - 212^ log 17 - 2f2|C - $ log fl\ 2 + 212K = 



j= - —+% + -A - ^x - m ■ xh + ny®c - Wiog n 

Ov 2 — 2 — — 

-0(-C + ^logO)S(-C + ^log!2)=0 

^-2^^-^|C| 2 + «C^log^ + ^|^log^| 2 + ^^K + itrxtrx-ix-xj = <> 



On C : 

9t 9 log Q 

-L-2n x + Cxl = Q , -^^+dxlogf2 + 2f2u; = 
d\ — o\ 

+ ^ tr X + dxtrx + 2fiwtrx + 12|x| 2 = 
|^ + trxC + Vtyvx - tytrx - ^° A g0) + O^logO-x + CxC- C x fi\og fi = 

+ fitrxtrx - 2fiwtrx + &trx - 2f2chVC - 212^ log 12 - 2f2|C + )W2| 2 + 2f2K = 

(/A 

9y fitrx 12try 9 log 12 

- -^x + -^x + —^x + (dxiog m - n( x ■ xh - n 7®c - n ^iog n 

-f2(C + ^log!2)§(C + ^logl2) +£^x = (2.51) 
^+7*" - 2f2^^--0|C| 2 - 12C^logl2 + -12|^logl2| 2 + -12f K + i trxtrx--x • xj =0 . 

Remark: The implementation of the initial conditions, namely the way of ob- 
taining initial data satisfying the constraint equations, \2.5(A \2.51\ with appropri- 
ate (Sobolev) regularity and asymptotic behaviour has been discussed in JCa-Nilf . 
Here we will have basically to repeat the same argument, but imposing the 
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analiticity, this makes this problem more complicated and how we solve it will be 
discussed later on. Next section, Section^ is devoted to find a local analytic so- 
lution for the system of equations \2.3T\ \2.38\ To do it we rewrite them in a more 
compact notation and show how, following G.F.D.Duff, JDuf and H.Friedrich, 
\Friy , we can apply the Cauchy-Kowalevski theorem to this characteristic case. 

3 The analytic solution of the characteristic prob- 
lem via the Cauchy-Kowalevski theorem. 

The method we use to obtain a real analytic solution of the characteristic 
problem defined by the system of equations 12.371 12.381 with initial data sat- 
isfying ecmations l2.501 l2".51l is a variant of the Cauchy-Kowalevski method. The 
adaptation of the Cauchy-Kowalevski theorem to characteristic problems has 
been developed by G.F.D.Duff, |Duj . for the linear case and, subsequently, by 
H.Friedrich, |Frlj . for the non linear problem. Friedrich result is suited to the 
present case, therefore we just recall the main lines of the proof, a straightfor- 
ward adaptation of his result. 

The system of equations l2.37ir2.38l can be written in a much more compact form 
in the following way: 

|^=F(V,W,^V,^W) 
dW 

^ = GV,W,fV (3.52) 
ov 

where V and W are vector functions valued in i? 18 and R 10 respectively, defined 

by 

V = {V s } = {^Q,v,t,^C,x} , W = {W*} = {X,w,u,x} 

s £ {1, 18} , t e {1,...,10} . (3.53) 

The initial data are assigned on the union of the null hypersurfaces Co and 
C_ . They have to be analytic functions satisfying the costraint equations 12.501 
12.511 this is possible as it is proved later on and has been proved in |Ca-Nil] for 
initial data belonging to a suitable Sobolev spaced Let us denote the initial 
data V = V (^,w a ), W = W (A,u a ) on C and C respectively^! 
The existence of a local real analytic solution of the system 13.521 with initial 
data Vo = Vo(i/,w"), Wo = Wo(A,w") is proved in the following theorem: 

Theorem 3.1 the system o f equations \3.52\ with initial data Vo = Vo(^,w a ), Wo 
Wo(A,w a ) admits a unique real analytic solution in a region (A, v) £ [0, A] x [0,F] 
whose size is determined by the initial data. 

24 The situation is somewhat simpler here with respect to Ca- Nil| due to the fact that we 
are considering a local problem and we do not have to worry about the asymptotic behaviour 
of the initial data, but only require that the initial data be analytic. 

25 The initial data defined here, Vo = Vo ( f, ui a ) , Wo = Wo(A,w°) satisfy the constraint 
equations 12.391 the remaining ones are automatically obtained using the Cauchy-Kowalevski 
method. 
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Proof: System 13.521 can be rewritten for the new unkown functions 
(V- V )(A,^,o; a ) and (W - W )(A, v, w a ) 



we denote again V and W. It is easy to show that it has the form: 



dV s _ wS , a dV s ' ^adW 1 ' 



(3.54) 




(3.55) 



with ui a G {8, 4>} and we sum over repeated indices. With these new V, W, the 
initial data are 



and the coefficients Fj/ a , F^/ a , f s , G*', a , g t depend, besides V and W, on the 
original data Vo and Wo and through them on the coordinates {x^}. 
Following Friedrich, |Frlj . the proof is basically made by two main steps. The 
first one provides a recursive mechanism to get all the derivatives in v, A, u a for 

V and W observing that, as initial data, we have all the derivatives in v, u a for 

V and in A, uj a for W. The remaining mixed derivatives are obtained through 
the equations 13.541 and 13.551 In the second and more delicate step we prove the 
convergence of the formal power series we have obtained. 

i) the recursive determination of the derivatives: From equation 13.551 
we control d v W and d^f q W for any q > 0. From equation 13.541 we control 
d\V and d\ft q V for any q > 0. Deriving equation 13.541 with respect to v we 
control d v d\fi 9 V, deriving equation l3.55l with respect to A we control d v dy^J q W . 
Deriving with respect to A equation 13.541 we control dfy/ q V and deriving with 
respect to v eauation l3.55l we control &lf q W . Iterating the procedure we obtain 
all the mixed derivatives. 

Remark: It should be clear that this procedure which allows to obtain formal 
power series for V and W both on Cq and on C_ satisfies also the constraint 
equations for V on <2 and W on Co as discussed in remark before Theorem 



ii) The convergence of the formal series 

The functions F s s / a , Fp a , G*', a , f s ,g l depend on the analytic initial data V and 
Wo and on the unknown functions V and W. More specifically, looking at the 
explicit expression of system 12.371 12.381 they are polynomials in V and W and 
can be written as 



V o = (0,...,0) , Wo 



(0,...,0) 



(3.56) 



EES 



s 



i.S'.a 





(3.57) 



a/3 




27 



with 



\a\ = , \/3\ < 2 for F , \a\ = , |/3| < 1 for F 

|a| = , < 2 for G (3.58) 

|or| < 1 , |/3| < 3 for / , |a| < 1 , |/3| < 3 for g . 

As we assumed that the initial data Vo,Wo are real analytic, the functions 
^(^) 5 ^;: /3 (^),G^(^),/^(^), 5 ^(^) are real analytic^ and can 
be written as power series in 



7 



G^W^E^,^ 7 (3.59) 



7 

( g 7 x . 



Due to the real analyticity we can assume that their convergence radius be 
> Rr^\ for a given R > and that in B R (0) C i? 4 the coefficients of the 
expansions 13.591 satisfy 

M 

\ F s';aBy\ ' \ F t'\aBj\ > \^s';aPy\ ' \fa/3y\ ' Iffa^l ^ • (3.60) 

Let us define the function H(8 fJ/ x' 1 ) = H(y), with 9\ > 1, 02 = O3 = 04 = 1: 

H{e ^) = _J"_ = = y ™ \f = y HyX , . (3 . 61) 

It follows from the previous definitions that 

l^s'';a/3 7 l ) \ F V°aPy\ > \ G a"ap<y\ > l/a/3 7 l > lfla/3 7 l < #7 ( 3 - 62 ) 

which means that the function H{9^x^) majorizes the functions F^' a a g(x), F^^x), 

G l'' a , a p( x ) ,fZp(x) ^ipix). We define now the functions F// a , #*; a , , f s , 5* 
as the power series |3~5T1 with the coefficient functions Fp a a p(x ,i )... substituted 
by the function H(9 fl x fJ ') times a matrix Pg! a a a-- with non negative real co- 
efficients and which defines the same polynomial r(V, W) = Yl a B PapW a V^ , 
PaB > , \a\ < 1 , |/3| < 3, for all indices s,s',t,a. Thereford^l 

F s s ; a = H{e^)yP^ ap W a V? = £T(<V)r(V,W)P s s /° 



26{ x /i} = { x i,a; 2 ,a: 3 ,x 4 } = {A,^w a }. 

27 As we are concerned about a local solution we consider a compact portion of Co UCq. 
28 The matrix elements P^, a , P^' a , Q ) a , p s , q l are all zero or one. 
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F s t ; a = H(6^)J2 P t'%P WaV ' 3 = H(6^)r(V,W)P t r 

a/3 

&f = H{6^) £ Q%pW a V? = H(0^)r(V, W)Q*f (3.63) 

a/3 

f s = H(6^) YttixfWV* = H(6^)r(V, W)p s 

a/3 

g l = H(8^)J2<lUW a VP = H(8^)r(V, W)q t . 

Let us consider the following system of equations: 

dV s - dV s ' ~ dW*' - - - 

BW t - f)V s ' 

^r = G - a ^ + ^ v ' w) ' (3 - 64) 

the following lemma holds: 

Lemma 3.1 The solution V, W of system \3. 6J\ with initial conditions Vo = 
Wo majorizing the initial conditions Vo = Wo = 0, is majorizing V and W, 
solution of system \3. 54\\ 3. 5 51 with initial conditions Vo = Wo = 0. 

Proof: The first step is to write formal power expansions for V and W; This 
is a formal expansion around a generic point of So = 5(Ao,fo)i redefining A 
and v as A — Ao, v — vq we consider it as an expansion around (A, v) — (0,0). 
The convergence of the formal expansion can be repeated for all the angular 
coordinates {uj a } of the points of So so that, finally, we obtain an analytic 
solution in a neighbourhood of So- Therefore, redefining also Lu a as (u a — cJq) 
we have: 

^) = E^ w T w a S ' W t (x) = J2™%a* j v k "i 1 "? (3-65) 

Observe that the initial conditions imply v^ ka = Wj 0a = and that, plugging 
the formal expansions in 13 . 541 13.551 we obtain recursively, as explained before, 
all the coefficients u| fca and Wj ka . We proceed in the same way expanding the 

solution V, W of system [3.641 
V'(x)=J2*'jka xil/ku ? u ? > W\x)=Y,™)kJv k uVu? (3-66) 

and we observe that, due to the definition of the functions F^' a , F s t ; a , G*', a , f s , g f , 
13.631 it follows that u| fea , W^ ka > and that the formal solution (V. W). 13.661 
majorizes the (formal) solution (V. W). 13.651 namely the following inequalities 
hold, for all index values, 

Vjka > \ v jka\ , wj% > Wjka\ • (3-67) 
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The proof is achieved if we can find a real analytic solution of the system 13.641 
with initial conditions such that the coefficients of its power expansion satisfy 
i'oka — i w t j 0a > 0; in this case, for an appropriate radius R, the expansions 
13. 651 describe a real analytic solution of system 13 . 5 5 1 with initial conditions Vo = 
W = 0. 

To find a solution of the svstem lB. 641 with initial conditions Vo = Wo such that 

«0fca>0,^0a >0 , (3.68) 
we make the following ansatz: 

V s (xn = V s (y = 9^) = V(y) 

W\x^) = W l (y = 9^) = W{y) . (3.69) 

Each equation for V s and W l of system 13.641 becomes in terms of V and W, 
denoting 

dy dy 

6{V' = H(y)q(V, W)(36V' + 20W' + c) 

W' = H(y)q{V,W)(36V' + 1) (3.70) 

where q(V, W) is a polynomial in V, W with positive coefficients; therefore 

(01 - 36Hq - 20 • 36{Hq) 2 )V' = Hq(20Hq + 1) 

W> = Hq + 36Hq + 



»i - 36Hq - 20 • 36{Hq) 2 
1 



[6xHq] (3.71) 



{Ox - 36Hq - 20 • 36(Hq) 2 ) 

and the two equations can be rewritten as 

(8 1 - 36Hq - 20 • 36{Hq) 2 )V' = Hq{20Hq + 1) 

{Ox - 36Hq - 20 • 36(Hq) 2 )W' = 9xHq . (3.72) 

Initial conditions for Vo = Wo such that 13.681 holds, are satisfied if we require 
as initial conditions for V and W 

V(0) = , W(0) = , (3.73) 

therefore we are reduced to prove that the system of equations 

(01 - 36Hq - 20 • 36{Hq) 2 )V' = Hq(20Hq + 1) 

(0! - 36Hq - 20 • 36(Hq) 2 )W' = O^q . (3.74) 

with initial conditions 13 . 731 has a real analytic solution in a neighborhood of the 
origin. Choosing 0i such that 

0i - 36# (0) g (0, 0) - 20 • 36(^(0)^(0, 0)) 2 > (3.75) 
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equations 13.741 can be written in the form 

V' = H(y,V,W) 

W' = J(y,V,W) (3.76) 

where H and J are real analytic functions having power series expansions at 
y = 0,V — 0,W — with non negative expansion coefficients. Therefore a 
real analytic solution of the system 13.641 exists and there is a neighborhood of 
So = S(\o,vo) where the real analytic solution of the system 13.551 does exists. 
This completes the analytic part of our result giving a concrete meaning to 
Theorem [2H1 

Remark: Recall that the local analytic solution whose proof has been sketched 
now is a solution for V — Vo, W — Wo- Then V, W are a local solution of the 
Einstein equations provided Vo and Wo are analytic functions satisfying the 
initial data constraints \2.5(A \2.51\ 

4 A common region of existence for the real an- 
alytic solutions of the non linear characteristic 
problem. 

We look for an "Alhinac type" result, |A-Mj . proving that, if some Sobolev 
norms, H s , of a real analytic solution are controlled up to a certain s0 then 
its "Cauchy-Kowalevski existence region" can be extended to a larger region 
whose size depends only on these Sobolev norms. This kind of results have been 
started by Lax [L], extended by Nirenberg |Nirj and proved, in the main lines, 
for the Burger equation by S.Klainerman and one of the authors (F.N.) |K1-Nil] . 
These results requires the hyperbolicity of the partial differential equations we 
are considering which is evident in the case of the Burger equation, but more 
hidden in the case of the Einstein equations F"l 

It is appropriate, before discussing our result for the Einstein equations, to 
look in a detailed way what has been obtained in the case of the Burger 
equation, |K1-Nil) . 



4.1 A summary of the analytic extended solution approach 
to the Burger equation in [Kl-Nilj. 

We give a complete survey of the various steps of this approach in the case of 
the Burger equation 

du du 

m +u o- x = °- (4 ' 77) 



29 For our problem they are explicitely defined in Subsection 1 1.21 

30 The hyperbolicity for the Einstein equations is explicit in the harmonic gauge, in a more 
general setting, for instance in the geometric gauge we are considering here, it expresses itself 
in the existence of a-priori estimates for some energy-type norms, see also Ch-Kl , Kl-Ni2 . 
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Step 1: We introduce a Banach space B a ^ p denned by the norm 



g \\D-f{;t)\\ L , n ^ n 



n=0 



where L 2 = L 2 (R n °) and, defining the multiindcx (3 = j3 2 , .. .,/?„„), 



n! 



i8;l/3|=n 



(4.78) 



(4.79) 



In the Burger equation case, no = 1, but we keep no > 1 as we are interested 
to extend this result to a more general case. We assume that the initial data 
of the hyperbolic equation we are considering are given by an analytic function 
which belongs to B atPo for a certain a > 2 and p > 0. 

Next lemma proves that if f(-,t) <E B a ^ p then, as a function of the "spatial" 
variables {x 1 }, f(-,t) is real analytic in R na . 



Lemma 4.1 Let f be a function belonging to the Banach space B ap with 

a > 



no 
2 



then f(-,t) is real analytic in B(x) p . 



Proof: If \\f(-,t)\\B atP < oo this implies that there exists a numerical sequence 
{h n }, depending on /, such that, for any n > 0, we have 

\\D n f(;t)\\^<^^, 
where {h n } is such that h n > and 

oo 

^ hn < 00 ■ 
n=0 

Therefore if / is a function e R n ° then by Sobolev Lemma, for \(3\ = n, 



(4.80) 



supj^/(x,t)i<co X! E n^/(-.*)iU a 

k=n P;\P\=k 



xeR n ° 



^ c o E /i; y^7 • 



n+l+[^l 

E 



kl h k 1 
n! fc Q 



(4.81) 



n! [(n + l)...(n+l + [^])] 

- V " ^ 

— n! n^"* 1 ^ — n! 1 
< c h n — — < c h n 



p n n a 



P" n 
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n 1 1 1 [, n °,1 

where h n = ^2 k = n 2 hk and, if p < 1, Co depends also on p. Therefore 

sup |D fl /(i,t)| < Co (n ,p)— ^— ^ ■ ( 482 ) 

x^R^o n( Q "l — 1) P n 

Given an arbitrary a; £ R n ° and y G B(x) p the formal power series 



(4.83) 



n=0 /3;|^|=n 

is convergent implying that in B(x) p , f(-,t) is a real analytic function. In fact 



ra=0£;|£|=n 



<c (n ,p)£ " ^ 



n=0/3;|/3|=n 

T 



]) ^ ft! ■ -ft" 



(4.84) 



As the choice of x was arbitrary and the estimate 14.821 is uniform in x it follows 
immediately that by unique continuation f(-,t) is real analytic in R n ° . In this 
case the neighborhoods B(x) p where the function can be expanded around x 
have also a common radius which is a stronger property. 
The previous estimate can be done in a better way we use later on, writing 



\D n f(;t)\\Li=Y,W Dn f(-^\\L 2 (Q(*i)) 



(4.85) 



where R n ° = Uj^ 1 Q(xj), where Q{xj) is a finite size square with "center" Xj. 
The following estimates hold 



|-D"/(-)*)IU 2 (0(sj)) 



< 



hn..iril 1 



(4.86) 



with h n j > such that 



h ^ = Z) Yl h ^ 



< CO . 



(4.87) 



n— j — 1 j — ln—0 

Repeating the previous argument we can conclude that if /(•, t) € B ap then 

n ! 



sup < Cl C/i nJ - 
leQfe) P 



(4.88) 

Next lemma says just the opposite of the previous one and its proof is immediate, 
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Lemma 4.2 Let f(t,x) be a real analytic function in a strip around the real 
axis of uniform height 2p, then for any T a < p f(t, •) belongs to the Banach 
space B a p i with appropriate p' < p and a. 

Step 2: The following theorem holds F*l 

Theorem 4.1 Let us consider ap.d.e. whose coefficients are real analytic func- 
tions, let vS°^ G B a ^ Po be the initial data defined on R n °; given p < po there 
exists a time T a (p, po, u^ ') and a solution u(t,x) £ C 1 ([0, T a ], B a p ) satisfying 
the initial conditions. 

This result follows more or less immediately recalling how Cauchy-Kowalevski 
theorem works; given (real) analytic initial data, using the p.d.e. equation 
it is possible to control all the derivatives in space and time of the (formal) 
series which defines the function u(t, x), analytic solution of our equation with 
initial data provided the series is convergent. Once this is proved the 
solution u(t, x) is real analytic with a convergence radius p, a priori depending 
on the point x S R n °] If belongs to -B QjPo this implies that the convergence 
radius of its power series is uniform in x. If the real analytic coefficients of the 
p.d.e. are defined, for instance, in the whole R n ° +1 with a common convergence 
radius, it follows that w(t, x) can be defined in a strip around the real axis of 
uniform width 2Ai. Therefore the solution exists for all \t\ < Ai. This does 
not prevent the possibility that choosing as initial data on the hypersurface 
with |f | = Ai — e the function u(t, x) it could be possible, again using the 
Cauchy-Kowalevski theorem, to extend the solution to a larger strip and then 
iterating this procedure forever or up to a moment when this has to end. The 
first case can happen, for instance, if the p.d.e. equation is linear as it follows 
that the strip width does not depend on t, while the second case takes place if 
in the iteration steps the width becomes smaller and smaller so that the strip 
where the analyticity is proved cannot be extended anymore. 
Nevertheless it has to be pointed out that in the statement of the theorem the 
time T a (p, p , w ') does not define the larger possible existence region, but only 
the region where the function u(t, x), thought as a function of the space variables 
only is real analytic, with convergence radius p. In other words T a finite does 
not prevent the possibility that u(t, x) can be extended as real analytic function 
in a larger region, but only that for t > T a the function u(t, ■) does not belong 
to Ba iP but, possibly, to a B a ^ p > with p' < p. To consider the (strip) largest 
possible analyticity region, possible infinite, we define 

T = sup{T Q >0\p(T a ,p Q ,u^) >0} . (4.89) 

If T < oo, it must happen that 

\im_p{T a , po) = (4.90) 

T a ->T 

31 Although this theorem is valid for a larger class of partial differential equations, we look 
at its result applied to an evolution equation of hyperbolic type. 
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and therefore in the limit t — > T the function lim. su(i, •) is not real analytic 
anymore in the space variables. 

The main conclusion from this discussion is that this approach does not provide 
a real control of the largest existence region. This can be a problem, for instance, 
if we are considering a sequence of real analytic solutions vS n ^ (t, x) associated 
to a sequence of real analytic initial data converging to a function belonging 
to a Sobolev space. In fact based on the previous theorem we cannot make a 
statement on the existence region of the limit solution as the can shrink to 
zero as n — > oo. 

In the case of the hyperbolic equation the situation is different. While in the 
general case, as discussed, the control of the time T is based, via the Cauchy- 
Kowalevski theorem, on the estimates of the coefficients of the p.d.e. equation 
and on the initial data u^°\ when the equation we are considering is hyperbolic 
one can have a better control of the solution existence time, namely one can 
prove that the analytic solution can be extended beyond T up to a time T 
which depends only on the suitable Sobolev H s norm of the initial data. In the 
following of this section we discuss in some detail the Klainerman, Nicolo result 
for the Burger equation, |K1-Nil) . 

We start considering the Burger equation solution u(t, •) <E C 1 ([0, T a ],B ai p) with 
initial data 6 B a ^ Po and with, as discussed before, 

P = p{Ta,Po,u {a) ) < Po ■ 
Energy conservation and Sobolev inequalities imply that 

\H;t)\\ H 2 < || u (°)|| ff2 expc / \\u(;s)\\ H 2ds (4.91) 

Jo 

and 

sup \\u(- 7 t)\\ H 2 < ||u^|| //2 expc / \\u(;s)\\ H *ds (4.92) 
se[o,T] Jo 

which implies 

sup [K,t)[[ff 2 < H^Hffsexp (cT sup \\u(;t)\\ H 2) . (4.93) 
te[o,r] te[o,T] 

Therefore we have the following lemma, 

Lemma 4.3 Let u(t,-) 6 C 1 ([0, T a ], Bq, iP ) be solution of the Burger equation 
with initial data £ B OiP0 , then there exists a time T > 0, depending only 
on \\u(°'\\ij2 such that for any T' < min{T,T a } 

sup \H-,t)\\ H i <c \\u^\\ m . (4.94) 
te[o,T'] 
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Inequality 14.941 is satisfied if T verifies, 

T ± — Irwj — ■ 

coc\\uW\\ H 2 

Remark: In the non linear hyperbolic equations it is possible that the apriori 
estimates do not hold beyond a certain time value due to the fact that the equa- 
tions do not have global solutions. In this case there will be a time value T* > T 
where some norm blows up. It could also happen that with more refined apriori 
estimates inequality like \4-94\ can be extended for all t values. In the case of the 
Burger equation we are in the first situation] 32 ]. 

Step 3: Next goal is to prove the following fundamental theorem: 

Theorem 4.2 Let w ' 6 B a pa be the initial data, then for t £ [0,T"] there 
exists an analytic solution, u{-,t) S -B Q , P " where p" < po and, moreover, p" 
does not depend on T a , but only on po, \\u^\\B a . Finally in the time interval 
[0,T'\ the following inequality holds with a constant c±, 

sup \\u(;t)\\ Bap „ <c 4 \\u^\\ B ^ m . (4.96) 
te[o,T>] 

Remark: The central part of this result is that p" appearing in the analytic 
solution 

does not depend on T a , 

P " = p"( P oM 0) \\B a ,J, 

while before, from the Cauchy-Kowalevski theorem, we had 

ueC 1 ([0,T o ],.B a , / ,) 

with p = p(T a , pq,u^). It is this fact that allows to extend the solution beyond 
the time T defined in \4-.89[ see the next theorem. 

Proof: The proof of the theorem is based on the following lemma, 

Lemma 4.4 Letu(t,-) £ C ([0,T a ], B ajP ) be the analytic solution of the Burger 
equation with initial data G B apo , with a > 2. Then, for any J, the 

following estimates holdW^ 

\\D J u(t)\\» < C^ J -^^\ , (4.97) 

J P'o 

where Cq is a constant satisfying 

co\\u {0) \\ Ba , PO <C (4.98) 



32 For the vacuum Einstein equations with small data we are in the second case, see ICh-Kl| . 
33 We use the simplified notation | \D J u(t)\\ L 2 = | \D J u(-, t)\ \ L 2 ■ moreover with D J we 



indicate the spatial derivatives, in the Burger case D J — 
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and Co is defined in the inequality \4-94\ 7 > 70, where 

70= sup |K-,t)|| H 2 (4.99) 
se[o,T'] 

and 7 — 70 depends on Co . 

Sketch of the proof of the lemma: The proof of Lemma l4~4l see (Kl-Ni), is 
inductive, therefore one starts proving that inequality 14.971 is true for J = 1, 20 
This is obtained from the apriori estimate 14.941 in fact, for J — 1,2 and t <T' 

2 

\\D J=1 ' 2 u{t)\\ L 2 < \\u(;t)\\ H 2 <c |K0)|| ff2 = C o^||W )|| i <4.100) 

Therefore it is enough to choose Cq such that 

\\D J=1 u(t)\\ L 2 <c |K0)|| H 2 < Ca- 
po 

\\D J = 2 u{t)\\ L2 < c \\u(0)\\ H 2 < C ^\ . (4.101) 

1 Po 

Both inequalities are satisfied if, with a constant c\ > co(2 q Pq), Co satisfies: 

C > Cl \\u^\\ H 2 . (4.102) 

We do not repeat here the inductive estimates, see pages 97,98 of [Kl-Nil] . but 
we recall that in the case of the Burger equation the structure of the estimate 
is 

\\D N u(t)\\ L 2 < e^+^To* I \D N u(0)\\ L2 + f dtF({D J<N u{t)}) (4.103) 



where F is a complicated expression, see eqs. (7.20),...,(7.23) of |K1-Nil) . which, 
nevertheless, depends only on the partial derivatives of u of order lower than TV 
and is estimated using the inductive assumption. To estimate ||I? Ar M(0)||i2 for 
an arbitrary N > 2 we can only use the fact that u(0, ■) = belongs to B ayPo 
and therefore satisfies 

||^(0)|| l2 <||u(°)|| Bq , po ^-L. (4.104) 

JV Po 

Therefore inequality [4T03] becomes, for TV > 2 and 7 > 570, 

N\ 1 rt 
JV Po Jo 



\\D N u{t)\\ L , < e ("^\\uW\\ B ^ Q ^± + I dtF{{D J < N u{t)}) (4.105) 



< e (^h*|| u (0)|| B _^^+ f dt F{{D^u{t)}) 

iy Po Jo 

< ^-^^ + f dtF{{D J < N u{t)}) (4.106) 

2 N a p" Jo 



In the General Relativity case the induction will start at a different value of J. 
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provided we choose Co satisfying, beside |4~1021 

Co > n^\\ Ba , P0 ■ (4.107) 

Using the inductive assumption the integration part in 14.1051 can be estimated, 
see |K1-Nilj . as 



dtF({D J < N u(t)}) < e (N-^Co^_ 1 



2C c 1 C2kik 2 
Pail ~ 7o) 



(4.108) 



where c±, c 2 , ki, k 2 are constants which do not depend on the solution norms. It 
follows that the lemma is proved if 

^dcafafal < j (4 log) 

Po(7-7o) _ 

Inequality 14. 1091 gives a lower bound on 7 — 70. Therefore Lemma 14.41 is proved 
provided Co, (7 — 70), ot satisfy 

Co > (d + 2)||u<°>|| B( , iPG , (7-70) > max(^ 2ClC ^ lfc2 C , 4 7o ) , a > 2 (4.110) 
Defining 

p' = poe-^' (4.111) 
it follows that the previous estimates for the J derivative become 

su Pte[0 , TI] \\D J u{t)\\L- < Co^^j • (4.112) 

j p 

and defining C3 such that 

c 3 |M 0) ||b q . po >C > Cl |hi( 0) || SQ , PO (4.113) 
we obtain, with p" < p 1 < po, 

00 / "\ N 1 

sup |Hi)lk p „<C £ [ P -\ < ca-^lltiWll < C4 ||^||b„. P0 ,(4.114) 
te[o,T'] ' JV=0 VP/ P ~P 



proving Theorem 14.21 . 

Step 4: We are left to prove that we can extend the solution beyond T a , this 
is the content of the following theorem: 

Theorem 4.3 Let u be an analytic solution of the Burger equation with initial 
data it(-,0) = it*- ' E -Bq.p with a > 2 i/ien the solution can be extended to an 
analytic solution E -B a ,p'" urat/i 

p'" < Pi = Poe~ 7T , (4.115) 

up to a time T which depends only on the Sobolev norm of the initial data 
\\h 2 > see \4-.95\ Moreover 

SUp \\u(;t)\\B apl „ <C 5 || U (0) || Bq , P0 . (4.116) 

te[o,T] ,p 
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Proof: Given the initial data G B a ^ Po we know that given T a there exists 
an analytic solution u such that for any t G [0, T a ] 



u{;t)eC\[0,T a ];B a 



pi 



with p = p(T a , po, u^ '). The previous theorem, Theorem 14.21 savs somethi ng 
much stronger, namely that there exists pi = poe~ lT such that for any T a < Tr 5 l 

«eC 1 ([0,T a ];fl a , w ) (4.117) 

where 7 and T depend only on ||u( '||b and ||u' ) ||jj2( R n ) respectively. Let 
us define T* in the following way: 

T* = sup{T a G [0,T] |for any t G [0, T a ] u{T a , •) G B a . Pl } . (4.118) 

As in the interval [0, T a ], u is an analytic solution in all the variables, it follows 
that taking the limit t — > T* also u(-,T*) is an analytic solution G B ap >>> with 
p'" < p\. In fact, as discussed before, T a does not define the largest region of 
the Cauchy Kowalevski solution, but the largest region where u(-,t) G B api . 
If T* — T we have proved our result. Let us, therefore, assume that T* < T, 
we can apply again the Cauchy-Kowalevski theorem defining the real analytic 
function ifW = u(-, T*) as the initial data on Yjt* and proving that there exists 
an interval [T*, T* + A] with T > T* + A such that, for t in this interval, there 
is a function, solution of the Burger equation with this initial data, 

v(t, •) G B a<p , 

with p > 00 

Observe that A depends only on p, pi and H'U^Hb,, p or, in other words, only 
p, po, T, \\u^\\b c , :P0 ■ Let us consider now the function w defined in the following 
way, for an arbitrary k > 1, 

w{t, ■) = u(t, •) G C fe ([0, T*],B a , p ,„) , for t G [0, T*\ 

w(t,-) = v(t r )eC k ([T*,T* + A],B a . p ) , for t G [T*, T* + A] . (4.119) 

It is clear that w(t, •) is C k in the time variable for t G [0, T* + A] and we can 
conclude that (p < p'"), 

ti)(V)6C k ([0,r + A] l B a ,;) 

and is a solution of the Burger equation in this time interval with initial data 
w(0,-) =u(°)(-). 

Let A be such that T* + A < T, we can apply again Theorem 14.21 to this 
function, prove that exactly the same inductive estimates holcf^l and conclude 
that 

w(t,-) G C fe ( [0, T* + A] , B a _ p n ) C C fc ([0,T* + A],B a , p ,„) . 



35 Recall that in the previous theorem p" < p' = poe~ lT and T' = min{T a ,T}, therefore 
if T a < T we can choose p" > pi = po e ' T - 
36 It follows simply applying Lemma 14. 21 
37 In this case T a = T* + A. 
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As T* + A > T*, it follows that T* cannot be the sup as defined in 14. 1181 unless 
T* = T. Therefore we conclude that with initial data u^ ' € B apg there exists 
a solution u in [0, T] such that 

H(t,-)eO*([o 1 nv) » 

which means that the size of the spacetime region where the Cauchy-Kowalevski 

solution can be extended, depends only on the H s=2 norm of 

Remarks: 

1 ) Observe that all the previous results are valid if the initial data belong to a 
different Banach space B ap i g with p' < po, the time T and the fact that the 
solution can be extended at least up to T does not change. 

2) The fact that Co is lower bounded by ||u^||s a p and not only on \\u^\\h2 
is not harmful and, on the other side, is what has to be expected. In fact, at the 
end of the whole proof, the various real analitic solutions u n (associated to the 
sequence of real analytic initial data {un ^}) will have the same existence time 
T, but their higher derivatives will have increasing bounds osn->oo. Therefore 
denoting Cq the constant Cq "associated" to iti°' ) it will follow that Cq"' — > oo 
as n — > oo . 

Step 5: the construction of a Sobolev solution. In this case we have 
a sequence of real analytic functions {un } converging in the H 2 norms to 
u (0) g fj2 j|- j g c i ear that due to this convergence the H 2 norms of all the 
functions are bounded by c||u^ ^ || h 2 ■ On the other side the L 2 norms of the 
higher tangential derivatives will in general diverge as n — ¥ oo, therefore we also 
have 

lim ||4 0) (-)Hb ( „, =oc (4.120) 

71— ^OO a ' p 

Then we have a sequence of real analytic solutions, {u n (t,x)} defined in the 
same interval of time [0,T] and the L 2 norms of the higher, > 2, derivatives 
satisfy 

\\D J u n {t)\\ L .<C { n) ±^- (4.121) 

3 Pin) 



and asC^ >c'||^ 0) (-)||s Q , p 



lim C^ n) = oo . (4.122) 



71— >0O 



This is exactly what we expect as the sequence of real analytic functions converge 
to a function which is not real analytic. The fact that this sequence has a radius 
p depending on n, follows immediately repeating the procedure to obtain 
the estimate 14.1211 and looking at the needed estimate 14.1091 which now reads 



2Co™^cic 2 fcifc 2 
Po(7 (n) - 7o) 



< 1 . (4.123) 



40 



As 7*™) has to increase it follows thalF^l 

P{n) < Pi"' = Poe^ ( " )T . (4.124) 

Therefore, as already said, although the radius of convergence of the sequence of 
the analytic solutions of the Burger equation, with initial data approximating 
Sobolev data, tends to shrink, nevertheless the time interval [0, T] where the 
analytic solutions do exist does not change. 

4.2 The extension of the previous result to the Einstein 
vacuum equations in the characteristic case 

The extension of the previous result to the Einstein vacuum equations in the 
characteristic case requires a detailed discussion, let us indicate all the differ- 
ences, with respect to the Burger case, we have to deal with. 

1) The initial data: The analogous of data on Eo for the Burger equation are 
the data given on C = Co UC . The discussion on how to obtain real analytic 
data satisfying the constraints on this null hypersurface requires again the use of 
a "Cauchy-Kowalevski-type" argument for the local part and a "Burger-type" 
argument for the global one. More precisely, let us consider first the initial 
data on Co, first we have to specify on So the tangential derivatives for the 
O and the O quantities, where only some constraints have to be fulfilled, see 
|Ca-Nilj . Then we use the transport equations, on (a portion of) Co, for the 
)P N O variables for any N. Next, once we have these estimates, we can use 
the transport equations for the ^7 O variables for any N. Finally we use the 
transport equations just as relations expressing derivatives with respect to v 
in terms of lower order v derivatives (in the case of the O quantities) or lower 
order v derivatives and tangential derivatives for the O quantities. Once all 
these estimates have been done we know that the O and O quantities satisfy 
the transport equations along (a portion of) Co which are exactly the "initial 
condition equations" of Theorem 12.11 equations 12.501 12.511 
Exactly the same has to be done on (a portion of) C inverting the role of the O 
and O quantities. This completes the local construction of the analytic initial 
data. 

Applying now a "Burger type argument" , we obtain estimates for all the deriva- 
tives in {v, uj b } of the O and O quantities proving that, being these quantities 
in a <B Q! p -type Banach space, they can be extended as real analytic functions 
on the whole Cq U C . The details of this construction are given in the following 
paper 

38 In the same way also in the construction of the initial data approximating Sobolev data 
we expect that po depends on n, po = Po 

39 The Burger argument developed before requires a bootstrap argument and, before, an 
apriori estimate. In the case of the initial data the mechanism is slightly different and although 
also in the initial data case, we need a recursive mechanism and a bootstrap argument, the 
fact that we can assign in a free way some connection coefficients makes the control of the 
initial data along the whole Co U C_ easier. 
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2) The local existence for real analitic solutions: Next step is to prove 
a local existence for our problem. This is done using the initial data on a local 
neighborhood of Sq and the version of the Cauchy-Kowalevski theorem adapted 
to this non linear characteristic problem as discussed in Section [3] This is the 
preliminary step needed to prove a larger existence region as discussed in Point 
3. 

3) The largest existence region: This is the central part of the proof. One 

has to show, as in the Burger case, that the existence region of an analytic 
solution depends only on some appropriate Sobolev norms of the analytic initial 
data. Arguing as in the previous discussion of the Burger equation this implies 
that we can have a sequence of real analytic data {\E , !°^} converging in some 
Sobolev norms to a "Sobolev" initial data whose corresponding solutions have 
a common region of analiticity. This would allow to construct H s solutions as 
a limit of a sequence of analytic solutions. 

The proof that the existence region of an analytic solution depends only on 
some appropriate Sobolev norms is obtained, as in the case of the Burger equa- 
tion, using a contradiction argument, namely we show that we can define a 
Banach space, £> Q . P , where the real analytic solution belongs (as function of 
space variables), when restricted to a certain region /C a , and prove that the 
solution remains in this Banach space even when the region is extended to a 
larger one, /C, whose size depends only on some Sobolev norms of the initial 
data. Of course, again, the choice of the Banach space must be such that if 
a function belongs to it and solve the ("hyperbolic") Einstein equations then 
it is real analytic and can be used as a real analytic initial data on the upper 
boundary of the largest JC a region to extend the solution to K.. Moreover this 
argument requires that the set of regions K, a is not empty which requires pre- 
liminary the local existence result discussed in 2). Therefore next problem is 
that of defining the Banach spaces B atP . 

4) The Banach spaces: To state precisely the analogy with the Burger equa- 
tion, the Banach spaces B a , P , has to be the analogous ones of B a>p , the region 
/C a the analogous of the region R x [0,T a ], the region K, of the region R x [0,T]. 
Let us recall the definition of the Banach space norm in the case of the Burger 
equation: 

IKV)IK P = E l|J " /( V )IU2 nV (4.125) 

n=0 

In the present case let W = ^>o + be a real analytic solution, where 

* = (V;W) = ( 7 ,n,i;,^w,C,x ; X,w,u,x) = ( 7 ,n,w,^0 ; X,w,0) 

*o = (rf,fl,v,il> ; X,w) , *i = (0;0), (4.126) 

in a region, {(A, ^)|(A, v) G [0, A a ] x [0, n a ]}, of the characteristic problem for the 
Einstein vacuum equations. Recalling that V denotes the part of the Cauchy- 
Kowalevski solution whose equations are in the incoming A-direction and W the 
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opposite one, we define 



\P(\ -)\\b., = E v zj 



n=0 



~ (sup pe{2 . 4} sup AeC ||/„(r, A;p)V"e>(V,-)||Lp(So)) 

o>,-)lk,,„ = E — - — i Lna p n > ( 4 - 127 ) 

where ||V™C(A, .)||lp(s )> I|V™O(iv)||lp(s ) are slightly symbolic expressions 
for 

\\V n O(v,-)\\ LP{So) := J2 \m? 9 0(v,-)\\ LHSo) 

p,q;p+q=n 

HV"0(A,-)|U,( So) := £ ll^0A,-)|U,( So ) > (4-128) 

ip(n) depends on which element of O is considered, f n (r, A) depends on which 
element of O is considered and has the structure /„(r, A) = r a ^~p |A| /3 '™', see 
for instance the analogous definitions for the first derivatives in |K1-Ni2] . Finally 

\\f(%-)\\i*(.s ) = (J g d^\f(\;v,uj a )\y , (4.129) 

where the measure is done with respect to a given metric tensor 70, analogous 
to the metric 70 on Sq, but with a different radial factor 

r (A,z/) = ro + ^-A) , (4.130) 

7o(A,,)(-,-) = ^^7o(v). (4.131) 

Finally if q is an /i-covariant tensor on Sq its pointwise norm \q\, defined with 
respect to the metric tensor %, isF"l 

\q\ 2 = qa ia2 ...aM lb2 ... b ^ a o lbi ro 2b2 ---r o sbs = \q\i ■ (4.132) 
The Banach space B a - yP is defined through the norm 

ll*(i)(A,^-)lk„ = HO(A, OIK, + Olls-p • (4-133) 

Remarks: 

?) T/ie norms 

sup sup ||r" + ^")-tv n O(A,.)|| LP(So) , sup sup \\f n (r, A;p)V n 0>, .)IU*W 
pe{2,4}^ec pe{2,4}AeC 



40 We use the metric 70 instead of the metric 70 as we want a "background" metric "near" 
to 7(A, v) for all the values of A, v. 
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play the role of the \\D n u(t, ■)\\L 2 (R n ) norm in the Burger case. The main differ- 
ence is that, in that case, D = ^ involves all the variables in R n (more specif- 
ically the only existing one, but in the case of R n with n > 1 all the derivatives 
{d x i} would be present). In the Einstein case the derivatives involved are those 
"tangent" to Cq and C_ repectively, namely {JI) 4 ,J7} on Cq, {tf> 3 ,y/} on C F^1 

ii) The previous norms are relative to the connection coefficients O and O while 
the Cauchy-Kowalevski solution ^ involves also the functions 7, f2, v, ip; X, w. 
Also for these quantities we have to define analogous norms; nevertheless once 
belongs to the Banach space B a;p we can easily prove that the analogous 
norms for these quantities are bounded. 

5) The main result: The core result we have proved and we present in the 
subsequent paper is that our solution belongs to a B a - P Banach space, with 
appropriate p and a, with (A, v) (coordinates of points) in a region K, whose size 
is determined only from some Sobolev norms |fl From this it follows immediately 
that in the same region the solution is real analyticlll Looking at the discussion 
of the similar result for the Burger equation we state, first of all, Lemma 14.51 
and Theorem 14.41 the analogous of Lemma 14.11 and Theorem 14.11 respectively. 

Lemma 4.5 Let 0(A, ■), Ofv, ■) be functions belonging to the Banach space B a - P 
with 

"31 



a > 



2 



then, for any x = {v,uj a }, 0(A, •) and, for any y — {A,u; }, Qfv,-), are real 
analytic in B(x) p , B{y) p respectively . 

Proof: The proof goes exactly as the proof of Lemma |4~T1 with the obvious mod- 
ifications and we do not repeat it here. Next result, the analogous of Theorem 
14. 1[ specifies the existence region of the analytic solution of the characteristic 
problem solved via the characteristic Cauchy-Kowalevski theorem as discussed 
in Section |3l 

Theorem 4.4 Let the "initial data" "J^ ) £ B a;Po then, given p < po, there 
exists a solution of the system of equations \2.37\ and ^(A, v, ui a ) 6 B a - P 

satisfying the initial conditions in a region, 

/C(A a ,n a ) = {(A, v) g [0,AJ x [0,n Q ]} , 

where A a ,II a depend on p,pQ,^>^ and are such that 

lim A Q (p,p ,* (0) ) = lim U a {p, Po , * (0) ) = . 

p->Po p->Po 



Proof: The proof goes exactly as the proof of Theorem 14.11 with the obvious 
modifications and we do not repeat it here. 



41 Applying p 4 ,^ 3 implies deriving with D4,D3 and projecting on the 5-tangent space. 
42 The region K is a region of the spacetime, but is completely defined once we give (A, v). 
43 In fact in a slightly larger region. 
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Theorem 14.41 tells us, exactly as in the Burger case, that if we define A, II such 
that 

/C(A,n) = sup{/C(A Q ,n a )|*(A,^,w a ) G B aiP with p > 0} (4.134) 

then 

lim_p(A a ,n Q ;p ,*(°)) = lim_p(A a ,n a ;p ,* (0) ) =0 . (4.135) 
A a ,— >a n a ^n 

Next step is to prove Lemma |4~61 and Theorem 14.51 which correspond to Lemma 
14.31 and Theorem 14.21 We first define the analogous of the H s norm for the 
Burger equation; we denote 

g=( 7 ,fi,X) , = (x,V,u) > 0=(x,2,w) 

and define, 

ll*(A,i/;Ollir.(So) = llsll^wM + ||0||jr.(s o )(A,i/) + \\Q\\a^s )M 

h\\ Hi (s ) = h\\ H;iSo )(\v) + \\n--\\ H}(So )(\v) + \\X\\ H s {So) (\u) . 
For the initial data we define analogous norms, 



v£C 

sup 



|* (0) ||UpEE Sup (||g||H.(S o )(0,I/) + l|O||fl-.(S )(0,I/) + \\O\\ Hl (s )M 

(l|g||^(s„)(A,0) + ||O|| 



Lemma 4.6 Let ^> — (g;0,0) = (7, f2, X; w, CiAi^x) ^ e a solution of the 
characteristic problem for the Einstein vacuum equations, in a region, 

£(A a ,II a ) EE {(A, I/) G [0,A a ] X [0,11a]} , 

wii/i initial data satisfying IH^^HLp < c> with p G {2,4}. Then there exists 
a region /C(A,IT) ; whose size depends only on the norm 1 1 1 \£>° 1 1 1 such that for 
s > 7, p € {2, 4} 7 there exists a constant Co such that 

sup ||M/(A^;-)||i«(5 ) <co|||* (0) ||U P (4.136) 
(x,i/)£/c(A',n') 

for any region K.(A',H') C min(K.(A, II), /C(A Q , II a )). 



11^1 



IIS 



[ff| = [J] / d^ |/ fe (r,A;p)y fc /(A;^a; a )r) 

\ i n ^ Sn / 
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Remark: Lemma \4-b\ states, basically, the apriori estimates for the Einstein 
equations. Differently from the simple Burger equation these estimates are much 
more involved as the hyperbolicity of the Einstein equations is more difficult to 
exploit. 

Next theorem, the analogous of Theorem 14.21 is the basic ingredient to prove 
our result. 

Theorem 4.5 Let = (g;0,<2) = (7, 0, v, ip, w, X; w,(,x,w,x) be a rea ^ 
analytic solution in the region, 

/c(A a ,n a ) = {(a, v) g [o,A a ] x [o,n a ]} , 

of the characteristic problem for the Einstein vacuum equations with analytic 

initial data on Cq U C_ q , belonging to the Banach space B a - Po . 

Then in any region K. C /C(A',IT) C /C(A a ,II a ) the following relation holds 

sup ||*i(A,i/;.)|| BaV , ^cJ^Hs^ (4.137) 

for some p" which depends only on po and on the | | |s Q . PQ norm, but does 
not depend on the region /C(A a ,II Q ) . 

Remark: Observe that p" and po are the analogous of the same quantities 
defined in Theorem \4-S\ in the Burger's equation case. Moreover (A a ,II a ) are 
the analogous ofT a , (A', II') are the analogous ofT' and (A, II) are the analogous 
of T and, exactly in the same way, they depend only on the apriori estimates 
which means that the size of the region /C(A,II) depends only on \\^^\\h s ' (S )- 
Lemma [4761 and Theorem l4.5l are the core of the result we want to obtain, exactly 
as the proof of Theorem 14.21 is based on Lemma 14. A\ the proof of Theorem 14.51 
is based on the following lemma, 

Lemma 4.7 Let = (g; 0,(2) = (7, v, -0, X; u, (, x, w, x) be a real ana- 
lytic solution in the region, 

K(A a ,n a ) = {(A, v) e [o,AJ x [o,n Q ]} , 

of the characteristic problem for the Einstein vacuum equations with analytic 
initial data on Cq U C_ , sufficiently small and belonging to the Banach space 
B a - Po , with a > 3. Then, for the generic connection coefficient we indicate with 
U, the following estimates hold for any J 

, j J\ P (J -2)5 p JT(\,v) 

\r 2+J - P f J U\ p , s (Kv) < C j^ -j (4.138) 

where 6 > Sq > 0, Cq is a constant satisfying 

Co > ci||*W|U M (4.139) 



45 The notation of inequality |4.138l is a bit symbolic, J is an integer except that in where 
it has to be cosidered a multiindex with \ J\ = J, moreover if with hi we denote a 0_ connection 

coefficient then the weight factor r 2+J v has to be modified, see equations 14. f 27l 
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and 



T(\,u) =T(u)+T(X) <7 (4.140) 
where, with an appropriate C > 0, 

T{v) = C ^A f r (A) = c(^2l. (4.141) 
vvq AA 

7 > 70 where 7 — 70 depends on Co and 

r (i/)+ro(A) < 7 o • (4.142) 

7o,<5o, To(^) + To(A) are relative to the initial data which are assumed to belong 
to B a - Pa and, moreover their angular covariant derivatives satisfy analogous 
bounds to \4-138\ with 7o,<5o, r (^) + r_o(A) instead of-f,6, T(is) +T(A). 

Remarks: 

a) The crucial fact in this lemma is that the estimates ^. 1 3 8\ do not depend on 
A a ,II a . Therefore the limit \J h .135\ is not anymore true and this basically allows 
to extend the region. 

b) Differently from Burger, when initial data are small a factor \v | + |A| in the 
exponential factor ( the analogous of T in the Burger case ) is not needed ( if we 
do not assume initial data "small" we expect nevertheless that the the \v\ + |A| 
has to be present). 

c) The goal of this lemma is to provide the appropriate estimates for the covari- 
ant derivatives of ^ to conclude that G -B Q ,p". Neverteheless to obtain this 
result we need to control the\ ■ | PjS norms not only for the angular derivatives, 
but also for the JZ) 4 derivatives for the V components and the the JZ) 3 derivatives 
for the W ones and also for the mixed derivatives. Nevertheless we first prove 
Lemma |^.7| and obtain the appropriate estimates for the angular derivatives, 
then using again the structure equations we control all the remaining mixed 
derivatives. 

Once Lemma [4761 and Theorem l4.5l are proved, next step is the proof of the final 
theorem, Theorem 14.61 , which is the analogous of Theorem 14.31 for the Burger 
equation and which we state here. Nevertheless it has be remarked that its proof 
is significantly different as we are considering now a characteristic problem. 

Theorem 4.6 Let * = (g;0,0) = (7, fi, V, ip, w, X; w,C, X>u,x) be a real 
analytic solution problem for the Einstein vacuum equations with analytic initial 
data on Co U C_q, belonging to the Banach space B a - Po . Then this solution can 
be extended to an analytic solution 6 B a]p in in the region /C(A,II) with 

p'" < Pl := poe- 7(|A|+|n|) (4.143) 

where 7 > depends on the B a - Po norms and A, II depend only on the Sobolev 
norms of the initial data. Moreover if the initial data are "sufficiently small" 
then the region where the analytic solution exists is unbounded ( in the v variable, 
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Proof: The first step to prove the theorem is the proof of the following lemma: 



Lemma 4.8 If the results of Theorem \4-4\ an d Lemma \4-6\ are true then it is 
possible to extend this solution to a new solution ^ real analytic in the region 
/C(A a + A,II a + A) and such that, for each A, v in this region, it belongs to the 
Banach space B a . p "' with, depending on the initial conditions, 

p'" < p e-T(l A l+l n l) . (4.144) 



Proof of Lemma l4.8t First we prove, using these propositions, that assuming 
as initial data those on Co U C and the real analytic solution on C(IT a ; [0, A a ]) 
and on C(A Q ; [0,IT a ]) we can extend the solution to a real analytic solution in 
/C(A Q , n a + A) U /C(A Q + A, n a ) and also on the diamond region 

/C(A Q + A, n a + A)//C(A a , n a + A) U /C(A a + A, n a ) . 

The proof of the existence of a real analytic solution in the two strips IC(A a , H a + 
A) and /C(A a + A, II a ) mimicks the analogous proof made in ;C a-Ni2 .. with the 
difference that, there, we were building an H s solution while here we prove 
the existence of a real analytic solution. The problem is nevertheless of the 
same type: as these stripes in one direction cannot have an uniformily bounded 
length we have to use a "sub-bootstrap" mechanism to prove the existence of 
the solution there. In the real analytic case the existence of the solution in a 
strip of a width which does not tend to zero is provided again by the "a priori 
estimates" which allow to show that the solution exists along the whole strip. 

Once Lemma 14.81 has been proved, the remaining steps to prove Theorem 14.61 
are basically identical to those used for the Burger equation, namely Theorem 
l4~3l Observe that instead of [0,T*] we have /C(A*,n*) and the definition of T* 

T* = sup{T a 6 [0,T] |for any t G [0,T o ] u(T a , •) G £? Q , Pl } . 

is substituted by 

IC(A*, TV) =sup{/C(A a , n a )c/C(A, n) | V (A, v) e /C(A a , n o ); *(A, v- •) g B a . pi } . 

Lemma T4.8I shows that, unless 

/C(A*,n*) = /C(A,n) , (4.145) 

the region K,(A*,H*) can be extended contradicting its definition. Therefore 
14.1451 is proved. 



5 Conclusions 

The results presented in this paper concern the existence region of the real 
analitic solutions of a class of characteristic problems for the vacuum Einstein 
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equations. As a byproduct to obtain it, we have discussed and developed a gauge 
associated to the double null cone foliation, first introduced in |K1-Ni2j . We 
consider this gauge very appropriate to this kind of problems and in particular 
we believe that to obtain the same global result proved here in a gauge like the 
harmonic gauge could be much more difficult, even if possible. 
The choice of the "double null cone gauge" has also the advantage, in our opin- 
ion, that it allows writing the Einstein equations as first order equations in a 
very natural and geometric wa^F^l because this foliation, differently from any 
spacelike hypersurfaces foliation, is more intrinsic due to the physical meaning 
of the null cones. Moreover in this approach the distinction, in the characteris- 
tic case, between those Einstein equations which can be considered "evolution 
equations" and those which have to be interpreted as "constraint equations" is 
completely clear. We also believe that our approach could be used to deal with 
the characteristic problem with "initial data" on a null outgoing cone hypersur- 
face, but this has not yet been worked. 

6 Appendix 

6.1 General aspects of the structure equations 
Proof of Proposition [2TT1 We have 

{tn A u°)(ep,e s ) = oj%{e s ) - uif(ep) (6.146) 

on the other side, recalling how d operates, 

d6 a (ep,e s ) = e f! (6 a (e s ))-e s (6'*(e ))-6 a ([ep,e s }) = -6 a {[ep,e s }) 

= -e a (T> eff e s ) + 6 a (D es ep) 

= -T} 5 6 a (e 7 ) + (e 7 ) = -T<* 05 + Tf p 

= -u%(ef})+u>%(e s ) = f? 7 Arf(ep,e 5 ) , (6.147) 

which proves the first structure equation. To prove the second structure equa- 
tion we write 

~Rj a /3 e s = R (ea,e/3)e 7 = D eQ (D e/3 e 7 ) - D e|9 (D eQ e 7 ) - D[ eQ e(3 ]e 7 

= D eQ (r^eA)-D e(j (r^e A )-D [eQie(3] e 7 (6.148) 
= e Q (r^ 7 )eA — e^(r^ 7 )eA + r^ 7 r^ A e e — r^ 7 r^ A e e — D[ eaje(3 ]e 7 

therefore 

R 7a/ 3 = e Q (r^ 7 ) - ep{T 5 ai ) + r^ 7 r^ A - r* 7 rjg A - 9 s (D[ ea , e ^e 7 ) 



46 This is certainly not new, see for instance |Ne-Pel| . [Ne-Pe2| . but we believe it is presented 
here in a more complete way. 
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On the other side 

(dw*+w* Aw^)(e Q ,e fl ) = dw$(e a , e ) + u 5 a {e a )u°{ep) - ui(ep)u^{e a ) 

= e a (uj s 7 (ep)) - ep(u}^(e a )) - ([e a , ep\) 
+ c4(e a K( e/3 ) - ui(ep)uj°(e a ) (6.149) 

and 

w*([e a ,e p ]) = ^(c^e CT ) = <f a ^{e a ) = c^T s aj 

= c°p6 s (D ea e 7 ) = 5 (D c ^e 7 ) = 6 s (D [eatefl] e 7 ) (6.150) 

and the thesis follows. 
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6.1.1 The general structure equations: du;^ = —uo s a A u° + 

P3X + trxx + 2uiX - f®Z + (2C - V - V)®£ = ~QL 

D 3 trx + ±{trx) 2 + 2utrx + |x| 2 - 2tyv{ - 2{- ( V + R - 2() = 

-curl£ + (2C + 3 - rj) A £ = 



D 4 irx + \trxtrx - 2wtrx + X ' X — 2c|lv?7 - 2|?7| 2 — 2£ • £ = 2p 
curlr/ + C A ^ - |"(x A x) = -ct 



p3X + \trx_X + \trxX ~ 2 ^X - f®V - V^>V ~ £®£ = 
D 3 fr X + \trxtr X ~ ^tr X + X ■ X ~ - 2|r/| 2 - 2£ • £ = 2p 

ci/rl?7 — \x A X + ^ A £ = a 



p4X + trxx + 2^x - f®(, - (2C + V + v)®£ = -a 

T> A tr X + \{tr X ) 2 + 2ujtr X + |xf - 2c}lv£ - 2£ • (5 + V + 2() = 

ci/rlC - (-2C + 77 — 2) A i=° 



cV rl X - C A x = 

ytrx - 4ivx + C ' X - Ctrx = -£ 



ci/rlx + C A x = — */? 

yirx - 4^ V X - C ' X + Ctrx = P 



p 4 | - J% -(r L - v )-x-4ui=-P 



p 3 ( + 2fu + x-(ri + 0-X-i+Mv-0-m = -P 
piC - 2fu - x(v -O + xC- Mv + O + 2cj£=-P 



ci/rlC- ^xAx = cr 

D 3 cj + D 4 u; - 4cjw + (7777 - £ • £ + • C — V ' = P 

2 Rdcab^dcAb + \trxtTX ~ \x ' X = ~P 
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6.2 The connection coefficients as functions on So X R x R 



To express the tensorial equations 12.311 12.321 as a system of p.d.e. equations 
we have to transform the covariant derivatives relative to the null directions 
into partial derivatives. This can be done expressing everything in a generic 
set of coordinates. We do it projecting all the connection coefficients on the 
two dimensional surface So = 5(0, 0) through the pullback associated to the 
diffeomorphism, 12.121 which specifies a coordinate set adapted to the double 
null foliation. 

Moreprecisely assuming the spacetime (M. , g) foliated by a double null folia- 
tion^ whose leaves we call, with a slight abuse of notation, outgoing or incoming 
cones and denote by C(A) and C_(v) respectively, the two dimensional surfaces 
S(X, v) = C(X) n Qiy) produce a foliation of each outgoing (incoming) cone, 
C(A), for instance, being foliated by the leaves {S(X, v)}. We recall that the 
vector fields 

N = 2fl 2 L , N = 2n 2 L, (6.151) 

are equivariant vector fields with respect to the leaves S(X. t/)F**l This means 
that the diffeomorphism generated by the vector field N, sends S(X, 0) onto 
$j,[5(A, 0)] = S(X, v), and $ A , the diffeomorphism generated by the equivariant 
vector field N_, sends S(0,v) onto $ A [5(0, v)] — S(X,v). We have previously 
defined, see 12.121 the diffeomorphism, 

*(A, v) : So 3 Po -> q = *(A, i/)(po) S S(X, v) (6.152) 

where *(A, v)(p ) = o $ A (po) = ®u(&\(po)) = <Mp) , V = $aM(6.153) 

From ^(A, v), we derive the pullback ^*(X : v) which sends the metric compo- 
nents and the various connection coefficients, tensors belonging to T*S(X, v) ® 

T*S(A, v) ® ® T*S*(A, v), to tensors belonging to T*S a ® T*S ® ® T*S 

and depending on the parameters A, v. We define: 

C(A,^;p ) = (**(A,i/)C)(po)=**(A J i/)(Co*(A J i/)(po)) 
= fi*X(i^AW)=ll$;((9) 

x(A,i/;po) = (**(A,i/) X )(po) = **(A,i/)(xo*(A,i/)(po)) 

= ®.*xKx($v o $ A (p )) = $aCx(<?) 

x(X,v; Po ) = (**(A,i/)x)(Po) = **(A,i/)(xo*(A,i/))(po) (6.154) 

= fj^X(** *a(Po)) = *A*Zx(ff) 

w(A, ^;p ) = (w o *(A, v))(p ) = (u o o $ A )(po) = w(t?) 
a)(A, ^;p ) = (^ o *(A, ^))(p ) = o o $ A )(p ) = cj((7) 

47 In |K1-Ni2l the double null foliation had to be "canonical" , see its definition in Chapter 
3. This is not needed here. 

48 Obviously they do not commute and their commutator is [JV, N] = — 4C 2 f (e\)^A ■ 
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where po 6 Sq = S(0, 0), is specified by its coordinates (9, </>). Therefore x, \-> C> 
Q, Q are covariant tensors defined on So and depending on the parameters A, v, 

x(A,i/;0,<£) , x(\v]6,4>) , l{\,v;Q,4>) , w(X,is;9,(j)) , u(X,is;9,(j)) . 

The Einstein equations will take the form of a system of first order partial dif- 
ferential equations, with respect to the variable {A, v, 0, (f>}, for the components 
of these covariant tensors and of the metric tensor. Defining on So an orthonor- 
mal basis {e^e^} with respect to a fixed metric 70 assigned on So, the explicit 
expression for £ is: 

C(A, v- 9, 0)(-) = Y, Ca(A, u; 9, <f>)6 A (-) = Ce(A, u; 0, 4>)d6{-) + &(A, u; 9, 0)#(-)(6.155) 

A 

where 

Ca(A, v; 0, 0) = C(A, v\ 0, 0)(e A ) (6.156) 

and exactly analogous expressions hold for x, x, w, w. Finally writing in .M the 
metric g(-, •) as 

g= \X\ 2 d\ 2 -2VL 2 {d\dv + dvd\)-X a {d\dLu a +duj a d\)+ lab dio a duj b , (6.157) 

we denote the pull back of its components, X,VL,^. In conclusion we have 
defined a diffcomorphism from (R 2 x So,g) to (.M,g) 

R 2 x So 3 p = (A, v, 9, 0) -> g = (A, i/)(po) (6.158) 

where po = </>) G So and 

g = |X| 2 dA 2 -2fi 2 (dA^ + d^dA)-X a (dAda; a + (iw a dA)+7 Qf) dw a da; fc . (6.159) 

The goal of this appendix is to rewrite equations 12.301 as a set of first order 
partial differential equations in the variables A, v, 9, <p. Following the discussion 
of Section 3 choosing an appropriate subset of equations 12.301 supplemented 
with some other equations for the metric components we obtain equations 12.311 
12.321 written (for the various components) in the {\,v,0,4>) coordinates. These 
are the Einstein equations written in the "double null foliation" gauge and 
Proposition 12.21 is. therefore, proved. 

6.3 The structure equations for the connection coefficients 
in a vacuum Einstein manifold in the {A, u, 8, <p} coor- 
dinates 

To write the previous equations in terms of the pulled back quantities we look 
first at those satisfied by £, corresponding to R(e4,e A ) = and R(e3,e A ) = , 

p4( + ^trxC + Cx - 4ivx + ^trx + P^logfl = 

Pzt + ^trxC + Cx + 4ivx - ^ftix - P 3 yiogft = . (6.160) 
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The derivative with respect to the parameter v of £ is 



(A,^;p ) = ($A < i ) :^C)(<z) ■ (6.161) 
av 



In fact 



i 

^(A, v; Po ) =1a Um £ C(*^-A ° IaCpo))" (K 0(<^ ° £a(po))] 

= iim («))- c(«)] = ®*:^o(*(a^)(po)) = mn^NOiq) ■ (6.162) 

A simple computation gives 

CCtfCXO = «[P4C(-) + (C-X)(-)] = fi[-tr X C(0 + (jlvx(-) - ytr X (0 - p 4 yiogfi(-)] 
which implies 

^(A,i/;M)(0 = (K*:"[-trxC + divx-ytr X -p 4 yiogfi])(A^;0,0)(-) . (6.163) 

Looking at the r.h.s. of 16.1631 it is easy to recognize that, denoting po a point 
of Sq the following relationships hold, whose proof is given later on: 

(*A#£fttrxC)(Po) = ^(A,^Po)trx(A,i/;p )C(A,^Po) 
CtA*^C-x)(Po) = ^(A,^;po)C(A,^;po)-x(A,^;po) 
(tI*Wvx)(Po) = n(A,i/;po)^x(A J i/;po) (6-164) 
(iI*^^ta:x)(Po) - 0(A,^;po)ytrx(A,^; Po ) 

$J*:np 4 yiogn(po) = ^fe^(A,i/;po)-nybgfi.x(A,i/;po) • 
In conclusion the pullback on 5o of the first equation of 16.1601 is 

3c ~ ~ ~ - — ~ (WiolTh 

+ fitrxC-^divx + ^ytr X + ^77^ f#logft-x = . (6.165) 

Of Of 

Equation 16.1651 has to be compared with the second equation of 12.301 The 
remaining equations along the outgoing cone are obtained in the same way. We 
recall that to write the pull back of the evolution equation along C(A) for x the 
following relation has been used @ 

Finally the remaining "outgoing equations" have the following expressions: 
dtr\ , ^trx~ , , Ai~i2. 



dv 



-trx + 2fiwtrx + fi|xr = 



9 The notation here can be misleading. With X'X w0 niean (x ■ x)ab = X ~1 cd Xdb> while 



with (x • x) we indicate (x • x) = X ai) 7 ac 7 M Xcd. 
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dtry ~~~ ~ ~ _•— -- _ ,„ 

+ Otrxtrx - 20u)trx - 2«4iv^ - 2ft|r)| 2 + 2flK = 

ov 2 — 2 ov — — — — — 

As expected from the fact that the vector fields N and N_ do not commute, 
the projection on 5*0 produces non equivalent expressions when applied to the 
equations along the incoming cones |£3 Let us consider again the equation for £. 
The derivative with respect to A of £ is: 

^(\,v ] Po) = *Z**u(CvQ(q) (6-166) 

where 

V = $*„AI (6.167) 
is the vector field generating the diffeomorphism ^o^o^ 1 . In fact, 
dC 1 

-± (A, u; po) = §X Km - [$tf>* u C($„ o §> x+h ( Po ))- ** C(*» ° *a (Po))] 

= Um ^ C**; 1 ^*! C(** ° *a+a(po))-C(^ ° *a(po))] (6-168) 

=t A $* hm i^v^^o^-ag)] =ti^(^c)(*(A^)(p )) =m:(£vO(q) ■ 

h-*Q h 

The following relation holds: 

{C v 0( q ) = (C K ()(q) - (£x()(q) (6-169) 
and £ satisfies the following equation: 



dX dX 



-fitrxC + ^#vx-^ytrx+yiogfi • x+£ x (- C xf l °g n = -(6.170) 



To express all the remaining equations on the incoming cones pulled back to Sq 
we use also the following equation, whose proof follows in the next subsection: 



®XK(ttp 3 f)(q) = ^(Po) + X\ P0 (f) . (0.171 ) 



Using also equation 16 . 1 7 1 1 we obtain 



dX 



^-2^-^|C?-^logft+^ = 



'Apart from the obvious substitution of A with v. 



55 



+ ntvxtrx - 2^wtrx + f x trx - 2Q$.vfj - 2ft|r?| 2 + 2QK = (6.172) 



dx ^ tr Xi fitrx- dlogfU 
9A ~ ~Y~ X + ~dX 



ritrxC + Oc^vx-oytrx- ^° g ^ + fayiogsi-x + - £^yiogO=o . 



+ ^7=trx + ^tr* + 2fiwtrx + ^|x| 2 = 



6.4 The structure equations in a vacuum Einstein space- 
time for the metric components in the {A, u, 6, 0} coor- 
dinates. 

To complete the set of p.d.e. equations representing the Einstein equations in 
this "double null foliation gauge" we have still to add some other equations. 
The previous equations are the analogue of the first order equations for the 
second fundamental form kij associated to the spacelikc foliation of the Einstein 
spacetime. We need the analogue of the equations for the three dimensional 
metric gij. Observe that in this foliation the metric written in the "adapted" 
coordinates is given in 16.1571 and the six quantity associated to the metric are 
Q, X a ,j a b. The corresponding quantities, pulledback to TSq, are 

7(A,^;p ) =K^7(^°£a(po)) =«a*Z7(«) 
X(X,iy; Po ) = 1aC*(^°$a(Po)) = ^* v X(q) 

Q(A,Iv;p ) = O($ 1/ o$ A (p )) = n(g) , 
and, proceeding as before, their partial derivatives with respect to v and A are 

0(A, v; 9, cf) = t A $*(£jv7)(-, ■)(?) = 2Ox0, '){\ "5 », <f>) (6.173) 

ov 



$1 
dX 



^(•, -)(A, v- 6, 0) =l^(ZV 7 )(-, ■)(?) = 2fix(-, ')(A, v\ 0, 4>)-(Cxl)(; •)(*> v, 0, , 



Beside these equations we need, in these coordinates, the equation connecting X 
to the connection coefficients and analogous equations for tt. From the definition 
of u> and w, in subsection 12. 3[ it follows 



u(\,u;po) = "T^ d 1 "^ (X,v;po) (6.174) 

u{X,v ]Pa ) =-i^^|^+y A lo g nJ(A,^p ) • 

The equation for X is derived from the commutation relation of N with N , 

[N,N] = -AQ 2 Z, (6.175) 
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It follows, in the {A, v, oj a } coordinates, 



dX 

— (\,is;p ) = -Q^^itfZ = 4n 2 Z(\,v;p ) , (6.176) 
av 



where Z{-) = Eag{i,2} C{^A)e A {-) ■ 

6.5 The equations 12.301 in the {\,v,9,(j)} coordinates. 

We collect all the equations 12.301 written now as pxl.e. equations for covariant 
tensors on So, depending on the parameters A, v f]] They are, indicating also 
the Ricci component to which they are associated, 



R(e 4 , e 4 ) =0 

trx + 20wtrx + il|x| 2 = 



atrx , fitrx— /, - o. 



dv 2 
R(e 4 ,eA) = : 

3? _ ^ - ~ „ (WlotrQ 

+ n tr X C - fidivx + nftrx + q - log fi • x = 

^AsR-(e J 4,es) = : 

c/fcrv — ____ - 

+ fitrxtrx - 2f2wtrx - 2f2div7) - 2ft|r)| 2 + 2fiK = 

R(eA,e B )=0 : 

<% fjtry » Otrv - <9 log f2 - - - - - - - - - 

^ - -^x + -^=x + —q^x - M{x ■ xh - nf®fi - O(2®2)=0 

R(e 3 ,e 4 ) = : 

^-2£iQu--n\(\ 2 +n(-fiogn+-n\fiogn\ 2 



R(e3,e 4 ) = 

<9cl> ^ - r,A~ - 3p, lIl9 a; a U,^,, a, 9 U/,-, 1~~ 1 

d\ x ~ 2 

<5 J 4sR(e J 4,e_ B ) = 



— +f^-2nQQ--n\c\ 2 -nc-fiogn+-n\f\ogn\ 2 +-n[ K+ 7 trxtrx--x-x 



+ fitrxtrx - 20u)trx + fxtr% - 2f2o/iv?7 - 2Q\fj\ 2 + 2QK = (6 

OA 



R(eA,e B ) = : 

9y Otrv » fitrv ~ <9 log Q » = ..... ~ ~, _ _ 
R(e 3 ,e j4 ) = : 



51 It is clear that this equations can be immediately written as equations for the various 
components of these tensors. In fact the orthonormal frame introduced in So does not depend 
on A, v. 
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|^ + n trxC + na/ivx - nftY K - + nfiogn-x + CjtC - £x?io g n=o . 

R(^e 3 ) = : 

^= + ^=trx + ^trx + 2^trx + ^|x| 2 = , 

where r\ = C + ^logfi and r\ = — £ + yiogi7 and R(-, •) is the traceless part of 
the Ricci tensor. To have a closed set of equations we have to add the equations 
at the level of the metric components. 16.1731 16.1741 and 16.1761 



?" 2 ^ = ° ■ 1 

OV OA 



2fix = , ^-2ftx + £x7 = (6.178) 



<91ogf2 - <91ogil si. ~ 

* + 2»£ = , -f +y y log» + 2»^= 

Ol/ OA 

^-4fi 2 Z = 0. 
a;/ 

6.6 The technical proofs. 
6.6.1 Proof of equations 16.1641 

($ A <i>*trxC)(Po) is the pullback of trx(<?)C(<7)- Therefore, denoting Yp a vector 
G TS'olp,, we have 

mK^xO(y P „) = m(^xO(q)(^^ x Y po ) (6.179) 

= (ntr X )(*(A,i/)(po))(tX$:0(Po)(5io) = ^(A,^;po)t?x(A^;p )C(A^;p )(r po ) . 
Proceeding in the same way we obtain 

m^KC-x){po) = n(A,i/;po)C(A,i/;po)-x(A,i/;po) • (6.180) 
To prove the third relation of 16.1641 

(*I$:n<Jfvx)(po) = n(A,i/;Po)^x(A,i';po) , 

observe that ckVx(-) is a one form defined on the spacetime; let Y po be an 
5*-tangent vector field at po, then 

mK^x)( P0 )(Y Po ) = divx(<?)($»^ A r P0 ) = <Jlvx(«)((*.„£. A 30 g ) 

= T^XKxWqX^*^^) = 7" 6 (9)(?j) fc (9)((<l',A/) ! ) e • (6.181) 



Immediately, with p = J> A (po), 

— i" c (Po) = $*,t* A F(p ) • (6.182) 
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Therefore 16.1811 becomes 



<9$f 



duj d 



duj c 



Y c (p 



Y c ( Po ) 



/'o 



so that, 



\dcu c 



IPo ■ 



(6.183) 



The remaining part of the proof is trivial and we do not report it herern To 
prove the last equation in 16. 1641 we observe that 

(OD 4 yiogfi)( eA ) = (C N flogn)(e A ) + (flogn)([N,e A })-flo g nCD N e A ) 
= {C N f\ogn-rtf\o g n-x){e A ). (6.184) 

Therefore 

^:np 4 f\ og n = $^*£^io g fi-$^*(nyio g n-x) 



dfiogn 

dv 



nfio g n- x 



(6.185) 



6.6.2 Proof of equation 16.1691 

We compute V — Q* U N_ in the following way: 

(®*uK)\ q (f) = N\ p (fo<f> v ) 



(6.186) 



where p = $ A (po) and q — <&v{p) — ^o^^po), /(•) is a scalar function defined 
on TM q and fo^JA is defined on TM P . Moreover p e S(X, 0) C C and we 



assume that on C r 3 l N 



J^. Therefore 



and from it 



K\p(f°*u) = ir 

ou 



($»*N)\ q 



dv 



df_ 

dxP 



9$, 



du 



d 

dxP 



(6.187) 



(6.188) 



As <& v is generated by the vector field N = J^, in the {A, v, to 11 } coordinates, a 
null geodesic along C(A) starting from the point p = (0, A, 9, <p) G S(X, 0) C C_ 



52 tt is enough to observe that concerning the indices a and b r y ab {g){^ a x)be{l) i s a scalar 
and therefore it can be rewritten in terms of the tilded quantities. 
53 As part of the initial conditions. 
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IS. 



$£(0, A, 9, 4>) = x±{v- 0, A, 6,4>) = v 

$£(0, A, 0> <£) = ^(f; 0, A, 9,4>) = A (6.189) 

*2(O,A,M)=z fl (i/;O,A,M) = 
$*(O,A,fl,^)=x*(i/;O,A,e,0)=^. 



This implies d$ v p /d\ = which, applied to 16.1881 gives 

d 



TM q 3V\ q = ($ v ,N)\ q = — 



= N\ 9 -X\ q (6.190) 



d\ 

and, therefore, 

(CvC)(q) = (%C)(<?) - 0CK)(«) ■ (6-191) 
6.6.3 Proof of equation 16.1701 

§(<?)(•) =1a<^ [(%C)(<?) - («)(?)] (•) =tA*I [(^ 3 c + oc-x-£xC)(g)] (•) 

= (*X«;n [(-trxC - d^vx + ytrx + JZ) 3 )P log - fi-^xC)] ) («)(•) 

where we have applied equation 16. 1601 satisfied by £ along the incoming cones. 
Let us compute $^$*np 3 yiogQ as we did in the previous case: 

(fiD 3 yiogO)(e A ) = %((yiogfi)(e^)) - ^log^D^) 

- (C K f log n)(c a) + (yiogfi)([Jv:,e A ]) - yio g r>(D i ve J 4) 

- (£jvyiogO - Oyiogr!-x)(eA) . (6.192) 

Observe that 

£xC = £*C (6.193) 

where 

X(p ) = X(q) . (6.194) 

Recalling eauation l6.1901 
dflogh 



dX 



C^flogfl - VlflogVL-x (6.195) 



54 Observe that X belongs to TSo, but has not to be confused with X on the points of So 
which is identically zero. 
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and, therefore the following relation holds: 
(tA$jn [(-trxC ~ 4^vx + ftrx + P 3 flog fi] ) (A, v; Po ) 

= ( - CitTxC - n^x + nfGx+ d ^ n + c^fhgn- nfi og n ■ x}(x,^ Po ) 

so that, finally, 

g _ ayiogn + ~ t ~ ~ + _ - + ^ log - _ _ + c _~ _ c ^ log n = 

6.6.4 Proof of equation 16.1711 

Let us consider |£ where f(po) — f(& v ° §L\(Po)) is a scalar function, 
df 1 

q^(pq) = lim ^ [/(** ° f/, ° ° *a(po) - /(^°$ A (Po))] 

= V|*„o* x (po)(/) = iV|($„o$ A (p ))(/) - ^l($„o$ A (p ))(/) ■ 

Therefore 

£*;W 3 /)(9) = = (OPa/)^ °*a(P0)) 

9/ df ~ ~ 

6.6.5 Proof of equations 16.1741 

ui(\,v;po) = °JbJ(Po) = ((D 4 logfi)o$„o$ A )(p ) 

' (Djvlognj^o^fpo)) 



2n($„o$ A ( Po )) 



' lim - [logfi($ ft o o $ A )(p ) - logr2($„ o t A )(p ) 



2fi(<I>„ o $ A (po)) ^->o ft 

Proceeding in the analogous way for w(A, u]Po) we obtain 

w(A,^;p ) = (w°iotA)W = -i((D 3 logr2)o$ iy o$ A )(p ) 
' -(£>jvlogO)(^o$ A (p )) 



2Q($ ! ,o$ A (p )) 

2»(^ 1 $ A ( Po )) I^logr!($, o $ A ( P0 )) + y x logfi($ w o *a(Po))] 
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' lim \ (log!J($„ o $ ft o 1 o o $ A )( Po ) - logl]($„ o$ A )( Po )) 



2f2($„ o $ A (po)) h ~*o h 



^logfi^o^^o)) = ___H_(A,i/;p )-y A logn(A 1 i/;po) 



6.6.6 Proof of equation 16.1761 

dX 1 

^(A,^;^) =$^ 1 $,; 1 lim-[$^ 1 X($, +h .o$ A (p ))-X($,o$ A ( Po ))] 

6.6.7 Proof of equations [2TMI 

Remark: Observe that 7(A, v){-, ■) is a symmetric covariant two tensor denned 
on Sq. We can consider it as a metric tensor, but the natural metric tensor on 
So is 7o( - j = 7(0, 0)(-, -)E On the other side, as discussed before, 7(A, v)(; •) 
is the pull back via <t*(A, v) of a metric tensor induced on S(X, v). The rela- 
tionship, reintroducing the tildas temporarily, is 16.1731 

|- 7 (A, !/;*„)(., •) = (¥*(A, i/)£* 7 )(¥(A, f)(po))(-, •) • 

Denoting {e^ } an orthonormal frame on So with respect to the metric 70 

it is clear that, in general, j(e^\ ^b) ^ $ab- The variables v cba (X, v, <J) 
are the (non covariant) derivatives with respect to the angular variables of 
7(A, v;po) — 7(A, v; oS) and do not have to be thought as the transport of partial 
derivatives of 7(A, v,Po) on S(A,v). The pullback is defined for the covariant 
tensors on S(X,v). This is true also for the covariant derivatives in the sense 

that denoted W a tensor on 5(A, v) it follows that fW = J/W and J7 is the co- 
variant derivative done with respect to the metric 7. Therefore if, for instance, 
W is a one form we have 

fW ab = f a W b = -^-W b - fl b W c (6.197) 



and T c ab = -j cd (v adb + v bad - v dab ) , (6.198) 



55 With respect to the metric 7o(-,-) = 7(0, 0)(-,-), wc define the area of So as | So 1 70 
1 

and its radius = |Sp| ^ n . As S p is diffeomorphic to S 2 , if the initial data on So arc 

"small", see for more details ICa-Nil . we can conclude that the metric 7o(',') is "near" to 
r (de 2 + sine 2 dcj> 2 ). 
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where v abc are the unknown functions of the first order system defined in 12.351 
Omitting the tildas, we have: 

g^ipa = g^f a log n - f a -^ log n = f a (-2fiw - y x b g n) 

= -2r»y7 Q ^ - 2Q^j a - (f a xy^ c - xy c ?p a . (6.199) 

d d 

~Q^ V cba = -Q^(d^)cba = d c —^f ba = -d c (C X l)ab + d c (2Q,x) ab 

= -{d c x d )d dlab - d x d clab + d c (d a x d ldb + d b x d lad ) + 2Qd c x ab + 2n{f c \o g n)x ab 

= -(d c X d )v dab - d x vcab + d c d a X b + d c d b X a + 2Vld c x ab + 2ttip c x ab 
= -(d c X d )v dab - d x vcab + d c w ab + d c w ba + 2ild c x ab + 2Vtij; c x ab 
To obtain the equation for w ab we write: 

dw cl b dd a X b d-/ bc X c 2 fp, d "fbc\ vc 

~5 — = —5 = d <* — 5 = - 4 <9a(^ Cb) + o a —— X 

dv dv dv \ dv J 

= -4d a (n 2 c b ) + d a (2n Xbc )x c 

= -8n 2 {f a log n)c b - 4n 2 d a c b + 2n{f a log n) Xbc x c + 2n(d aXbc )x c 

= -8QVaC& - ^daCb + 2^ aX bcX C + 2n{d aX bc)X C . 

We collect the three equations obtained, 

^-4 a = -2no a u - 2nu»p a - (f a xy^ c - xyjj a 

d 

^w c fca = -(d c X )v dab - d x v ca b + d c w ab + d c w ba + 2Qd c x_ ab + 2£l-ip c x ab 

9w ab _ Qo2 „, , ,o2 



dv 



= -80 2 vUb - 4n 2 d a c b + 2n^ aXbc x c + 2n(d aX bc)x c . (6.200) 



6.7 The null components of the Riemann tensor on So x 
R x R and their Bianchi equations 

In our notation the ten independent components of the (conformal part of 
the) Riemann tensor are denoted a, /3, p, a, /3, a, respectively two covariant two- 
tensors, two covariant vectors and two scalar functions, tangential at each point 
to S(X, v), the two dimensional surface containing it. All these components writ- 
ten in terms of the connection coefficients and their derivatives are, see |Ch-Kl] . 
|K1-Ni2j : 

a = -piX - trxx + (D 4 log fi)x 
/3 = ftrx - divx - C ' X + C tr X 

p = -K- -trxtrx + -\ ■ X 

a = cirlC- ixAx (6.201) 
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§_ = - ftrx + 4^ v x - C • X + C tr X 

QL = -p3X - trxx + (D 3 logft)x 

We will consider these quantities pulled back on the two dimensional surface 
So, as we did before for the connection coefficients. Therefore we define 

a{\v;p Q ) = (*"(A,i/)a)(po) = ° £a(po)) = *a$X«) 

^(A^; Po ) = (**(A,i/) j 8)(po) = *X#:^o$ A ( Po )) = *J*; ) 8( g ) 

£(A,»/;po) = (tf*(A,i/)£)(po) =ll$m°fA(l'o)) =W 

a(A,^;p ) = v)a)(po) = lS$*a($„ *a(po)) = 

p(A,i/;po) = (po*(A,i/))(po) = (p°$i,°$a)(Po) =p(g) 

er(A, i/;po) = (o- ° *(A, i/))(po) = (a o o t A )(p Q ) = <r(g) (6.202) 

where po € So and is specified by its coordinates (9, cj>). In conclusion we have 
the following covariant tensors defined on So, 

a(\,v;0,4>) , a(X,v;6,4>) , j3(\,i>;6,<j>) , j3(\,v; 6, <f>) , p(\,v;6,(/>) , a(\,v;6,</> 



Lemma 6.1 The explicit expressions for the pullback of the various null Rie- 
mann components are 

a = + x 7 - 2wx 

Q av 

P = -jf tr x - 4 iv x - C- X + 2^ trx 

p = -K- ^trxtrx + ^X'X 

o=cTrll-\x^X (6.203) 

§_ = ~ftrx+ (fivx-(-X+ \c tr X 
1 / dx 



Proof: It follows immediately from the explicit expressions of the various Rie- 
mann null components, 16.2011 and the pullback of the various connection coef- 
ficients and their derivatives. 5 3 

The Bianchi equations take the form of a system of transport equations with 
respect to the variables A and v for these covariant tensors. Defining on 5 an 
orthonormal basis {eg, e^}, each of them can be written in the following way: 

&(■) = &ab(X, v, e, <P)8 a ®0b(; •) (6.204) 

A 
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where oiab{\ v\ 0, 4>) = a(X, u; 9, <p)(e Al e B ) ■ 

We start writing the derivative with respect to v of a obtaining: 

^(A,i/;p )=^$;(£j V o)(g) (6.205) 

The proof goes over exactly as the one for §£(A, v]Po) and we do not repeat it. 
A simple computation gives 



(C N a)(-) = np4tt(-) + n(a- X + X-a)(-) = Q [p 4 a + tr X a + (a-xfr] (•) 



+ (a • x)7 - + (4^a - 3(xp - *X°) + (C - 4??)§^) 



(6.206) 



which implies 

g(A,M,40(-) = (± 'I' " 



-^a + (a • x)7 ~ f®§_ + (4^a - 3(*p - *X^) 



+ (C-4r?)§^) 



(A,i/;M)(0 



(6.207) 



Applying, as before, the pull-back on the right hand side, the final expression 
of the Bianchi equations, for this component, is 



da 
1h> 



(A, i/; M)(0 = « 



+ (C 



tr X _ , 

—a + (a • XH 



^4wa - 3(xp - 



) 



(A,i/;M)(0 . 



(6.208) 



It is easy now to project the remaining Bianchi equations on So. 

dp 



+*fa-3(rip- V)])] 



dp 



da 
dv 



-X-« + C-/3-2 > 71ogfi./3 



= fi(--tr X/ 5 + <jM- 

= $i$:(r»D 4( x) = $*<i>: 



0( __ trX a + div*/3 



-X-a + C-/3-2y71ogfi-^ 



"I --trxa + ckV*/?- 



-X-*a + C-*/3-2yio g r!. 



-X-*a + C-^-2yiogQ-^ 



n (-2 tr X+44va-2cj/3+x-/3+(C + yiog^)-a 



fi( -2trx/3+4iva-2o;/3+[x^+(C+yiogO)- Q 



(6.209) 
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The Bianchi equations which are transport equations along the incoming cones 
can be obtained exactly in the same way as done for the connection coefficients. 
In this case we have to use the relation, see 16.1691 for a generic Riemann null 
component w, 

^(A, v- 9, </>){■) = mt [(£Nw)(q) - (Cxw)(q)] . (6.210) 



d§_ 
dX 
dp 

dx 

da 
dX 
dp 
dX 

da 
dX 



nl--trxl3 + x-f3_-4va- 2ujfi + {-2( + rj) ■& 



1 



(6.211) 



(a-xh + f®P + (4^5 - 3{xp + Iff) + (C + 4?7)®/3) 
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